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CONVERSION  FACTORS,  NON-SI  TO  SI  (METRIC) 
UNITS  OF  MEASUREMENT 


Non-SI  units  of  measurement  used  in  this  report  can  be  converted  to  SI 
(metric)  units  as  follows: 


_ Multiply 

cubic  feet 
degrees  (angle) 


inches 

pounds  (mass) 

pounds  (mass) 
per  cubic  foot 

square  feet 


0.02831685 

0.01745329 

0.3048 


0.4535924 

16.01846 


0.09290304 


To  Obtain 

cubic  metres 

radians 

metres 

centimetres 

kilograms 

kilograms  per 
cubic  metre 

square  metres 


STABLE  RIPRAP  SIZE  FOR  OPEN  CHANNEL  FLOWS 


CHAPTER  1 
INTRODUCTION 

The  transport  of  water  through  natural  and  man-made  open  channels 
carries  the  possibility  of  scour  if  the  channel  boundaries  are  erodible. 
While  many  different  methods  have  been  used  to  protect  channel  bound¬ 
aries,  riprap  revetment  continues  to  be  one  of  the  most  widely  used 
methods.  Riprap  is  long-lasting,  flexible,  easily  placed  and  repaired, 
and  natural  in  appearance.  However,  in  some  locations  riprap  is  not 
readily  available  or  the  available  stone  is  too  small  for  riprap.  In 
other  locations,  a  limited  number  of  available  gradations,  rather  than 
design  guidance,  determines  the  size  used.  Transportation  costs  for 
riprap  from  quarry  to  jobsite  are  often  greater  than  the  cost  of  the 
rock  alone.  In  spite  of  these  limitations,  the  large  amount  of  riprap 
used  requires  guidance  to  ensure  optimum  design. 

Determining  riprap  size  is  one  of  the  most  important  factors  in 
defining  the  optimum  riprap  gradation.  Existing  riprap  sizing  methods 
have  limitations  which  include  the  following: 

1.  Many  existing  riprap  sizing  methods  have  evolved  from  sediment 
transport  concepts  which  use  shear  stress  to  define  particle 
stability.  Critical  shear  stress  for  a  given  riprap  size  is 
determined  by  the  well-known  Shields  coefficients.  Most  sedi¬ 
ment  transport  and  riprap  sizing  techniques  use  a  constant 
Shields  coefficient  for  rough  turbulent  flow.  Existing  riprap 
design  techniques  also  use  logarithmic  velocity  laws  to  relate 


2 


velocity  to  shear  stress.  However,  several  investigators  have 
found  the  Shields  coefficient  to  vary  at  high  relative  rough¬ 
ness  while  others  have  found  the  logarithmic  velocity  laws  to 
be  affected  by  high  relative  roughness.  Since  most  riprap 
stability  problems  involve  high  relative  roughness,  many  of  the 
existing  riprap  sizing  methods  may  not  be  applicable. 

2.  Existing  riprap  sizing  methods  that  use  shear  stress  have  an 

additional  liability.  As  stated  by  Neill  and  Hey  (1982), 

Researchers  tend  to  favor  shear  stress  criteria  for 
stability  and  bed  movement.  From  a  practical  engineer¬ 
ing  viewpoint,  local  shear  stresses  are  difficult  to 
measure  and  to  conceptualize,  compared  to  velocities. 
Researchers  might  pay  more  attention  to  expressing 
results  in  velocity  terms  for  practical  applications. 

3.  Existing  riprap  sizing  methods  also  lack  variation  relative  to 
cue  effects  of  riprap  gradation,  thickness,  and  shape. 

4.  The  analytical  techniques  used  to  determine  the  decrease  in 
stability  that  results  from  placing  riprap  on  a  channel  side 
slope  need  to  be  tested  against  experimental  data. 

Considering  these  limitations  of  existing  riprap  sizing  methods, 
the  objectives  of  this  study  are  as  follows: 

1.  Evaluate  the  applicability  of  existing  riprap  sizing  methods 
that  use  a  constant  Shields  coefficient  or  the  logarithmic 
velocity  laws. 

2.  Develop  a  riprap  method  based  on  velocity.  Determine  which 
velocity  (bottom,  average,  surface,  or  maximum)  to  use  in  the 
riprap  sizing  method. 

3.  Incorporate  riprap  gradation,  thickness,  and  shape  variation 
into  riprap  sizing  method. 


4.  Evaluate  side  slope  effects  on  riprap  stability  and  incorporate 
into  riprap  sizing  method  for  straight  and  curved  channels. 

A  series  of  flume  tests  were  used  to  accomplish  these  objectives  by 
studying  the  stability  and  resistance  to  flow  of  riprap  having  grada¬ 
tion,  thickness,  and  shape  similar  to  that  used  for  scour  protection  in 
open  channels.  Results  are  limited  to  channels  with  slopes  less  than 
2  percent,  and  the  ratio  of  flow  depth  to  average  riprap  size  must  be 
greater  than  4.  Riprap  sizing  for  placement  in  highly  turbulent  flow 
downstream  of  hydraulic  structures  or  for  placement  on  embankments 
subject  to  overtopping  flows  is  not  covered  in  this  study. 

The  following  chapters  present  first  a  review  of  existing  litera¬ 
ture  relative  to  these  four  objectives.  Next,  the  experimental  investi¬ 
gation  is  explained,  and  then  the  analysis  and  results  to  achieve  each 
of  these  four  objectives  are  presented.  Finally  the  conclusions  from 
the  study  and  recommendations  for  further  studies  are  presented. 


CHAPTER  2 


LITERATURE  REVIEW 

In  the  study  of  open  channel  riprap  stability,  many  investigations 
have  been  conducted  that  are  applicable  to  this  engineering  problem. 
This  review  of  existing  information  focuses  on  four  different  topics 
which  correspond  to  the  four  objectives  in  the  Introduction.  First, 
studies  concerning  the  effects  of  relative  roughness  on  Shields  coef¬ 
ficient  and  logarithmic  velocity  equations  will  be  reviewed  to  see  if 
existing  sizing  techniques  are  valid.  Second,  the  literature  will  be 
searched  for  existing  riprap  sizing  methods  based  on  velocity.  Third, 
previous  studies  will  be  reviewed  to  determine  the  present  knowledge 
regarding  the  effects  of  thickness,  shape,  and  gradation  on  particle 
stability.  Fourth,  existing  concepts  of  side  slope  particle  stability 
will  be  reviewed  and  summarized. 

2.1  APPLICABILITY  OF  EXISTING  RIPRAP 
SIZING  METHODS  USING  A  CONSTANT 
SHIELDS  COEFFICIENT  OR  THE 
LOGARITHMIC  VELOCITY  LAWS 

One  of  the  most  common  methods  for  evaluating  riprap  stability  is 
the  critical  shear  stress  method  (also  called  tractive  force).  The 
shear  stress  stability  concept  was  used  by  Dubuat  (1786)  but  did  not 
become  popular  until  Schoklitsch  (1914).  Lane  (1953)  used  the  tractive 
force  method  for  stable  canal  design  in  noncohesive  material.  Anderson 
Paintal,  and  Davenport  (1968)  developed  the  tractive  force  approach  int 
a  riprap  design  method  which  includes  the  effects  of  side  slopes  and 


channel  bends.  The  work  of  Anderson,  Paintal,  and  Davenport  is  used  as 

the  basis  for  riprap  design  by  the  US  Department  of  Transportation 

(1975).  The  Office,  Chief  of  Engineers  (OCE) ,  (1970  and  1971)  riprap 

design  guidance  is  based  on  the  tractive  force  approach.  Li  et  al. 

(1976)  and  Stevens  and  Simons  (1971)  developed  tractive  force  methods 

which  incorporate  probability  and  safety  factors  into  the  design  method. 

The  shear  stress  exerted  on  the  boundary  in  uniform  flow  is 

T  =  Y  DS  2.1 

'w 

where 

T  =  tractive  force  imposed  by  flowing  water* 

Yw  =  specific  weight  of  water 
D  =  flow  depth 
S  *  energy  slope 
or  using  hydraulic  radius 

T  =  y  RS  2.2 

'w 

where  R  is  the  hydraulic  radius. 

The  imposed  force  calculated  from  either  Equation  2.1  or  2.2  is 
equated  to  the  ability  of  the  particle  to  resist  movement  or  the  criti¬ 
cal  shear  stress.  Using  the  analysis  of  Carter,  Carlson,  and  Lane 
(1953),  which  is  an  equilibrium  force  analysis,  yields 

T  =  C.  (y  ~Y/d  tan  $  2.3 

c  1  \  '  s  'w  / 

where 

Tc  =  critical  tractive  force  for  given  particle  size  on  bottom 

*For  convenience,  symbols  and  unusual  abbreviations  are  listed  and 
defined  in  the  Nctation  (Appendix  A) . 


Cj  =  coefficient 


=  specific  weight  of  stone 
d  =  particle  size 
<5>  =  angle  of  repose 

Formulations  of  the  shear  relations  from  dimensional  analysis 
depend  on  which  parameters  are  considered  significant.  Vanoni  (1977) 
uses  the  parameters  T£  ,  “  Yw  *  d  »  the  fluid  density  p  ,  and 

viscosity  v  ,  to  define  incipient  motion.  This  results  in  the  same 
form  derived  by  Shields  (1936)  or 

T 

(Ys  -  Y»)d 

where 

U*  =  shear  velocity 
g  =  universal  gravitational  constant 

u*d 

For  rough  turbulent  flow  (particle  Reynolds  number  — —  >  400),  the 
right  side  is  often  assumed  constant  and  called  the  Shields  number  or 
Shields  coefficient,  herein  denoted  as  .  Most  of  the  stability 
investigations  concerning  Shields  coefficient  have  been  related  to 
sediment  transport.  According  to  Graf  (1971),  the  definition  of  the 
critical  Shields  coefficient  has  been  subject  to  the  interpretation  of 
the  researcher.  The  riprap  design  procedures  by  OCE  (1970)  and 
Anderson,  Paintal,  and  Davenport  (1968)  use  a  constant  Shields 
coefficient  for  safe  design. 

The  use  of  a  constant  Shields  coefficient  has  been  questioned  by 
Meyer-Peter  and  Muller  (1948),  Yalin  (1965),  Barr  and  Herbertson  (1966), 
Blench  (1966),  Neill  (1967  and  1968),  Bogardi  (1968),  Ashida  and  Bayazit 


(1973),  Bathurst,  Graf,  and  Cao  (1982),  Daido  (1983),  and  Bettess 
(1984),  who  propose  that  Shields  coefficients  should  vary  with  relative 
roughness.  Bathurst,  Graf,  and  Cao  (1982)  and  Bettess  (1984)  have  found 
this  variation  with  relative  roughness  to  be  limited  to  high  relative 
roughness  below  which  Shields  coefficient  becomes  constant.  Meyer-Peter 
and  Muller  (1948)  found  that  the  limiting  shear  stress  is  proportional 
to  particle  diameter  and  relative  roughness  and  proposed  an  equation 


An  explanation  for  a  changing  Shields  coefficient  with  relative  rough¬ 
ness  has  been  offered  by  Escoffier  (1968).  At  high  relative  roughness 
(low  depth/d^^),  turbulence  generated  at  the  boundary  is  hindered  by  the 
presence  of  the  free  surface.  Consequently  the  fluctuations  in  velocity 
are  decreased.  At  low  relative  roughness  (large  depth/d^),  the 
boundary-generated  turbulence  is  not  hindered  by  the  free  surface  and 
fluctuations  in  velocity  are  not  reduced.  Since  the  magnitude  of  turbu¬ 
lent  fluctuations  is  critical  for  riprap  stability,  this  provides  an 
explanation  for  the  variation  of  Shields  coefficient  with  relative 
roughness.  Chen  and  Roberson  (1974)  and  Bayazit  (1976)  found  that  mea¬ 
sured  turbulence  intensity  decreased  with  increase  of  relative  roughness 
in  the  region  near  the  wall.  Bayazit  (1982)  proposed  that  this  "can  be 
explained  by  the  fact  that  a  substantial  part  of  the  energy  of  the  mean 
flow  is  converted  into  turbulence  in  the  separation  zones  between  the 
roughness  elements  in  the  case  of  large  scale  roughness."  Gessler 
(1971)  stated  that  relative  roughness  does  not  influence  Shields  coef¬ 
ficient  because  incipient  conditions  depend  only  on  conditions  at  the 
bed  and  not  on  the  boundary  layer  thickness  (or  depth  in  open  channels). 


' j>  v  v  wavjWCj t/w. v\  v.  c-  w'*.-  v  /vuv -jv ..v_w»i _v  Awmv.v." « v*  w  ■  ■' 


Some  of  the  existing  riprap  procedures  (OCE  1970  and  Li  et  al. 

1976)  use  the  logarithmic  velocity  laws  to  determine  the  relation  between 
velocity  and  shear  stress  on  the  boundary.  The  universal  velocity  dis¬ 


tribution  law  for  rough  surfaces  is 


y  _  2.3 


30(y  +  yQ) 


where 

=  local  velocity  at  distance  y 
<  =  von  Karman  coefficient 
y  =  distance  above  origin 

yQ  =  distance  below  top  of  roughness  element  to  origin  of  profile 
Ks  =  equivalent  sand  grain  roughness 

Equation  2.6  is  integrated  over  the  depth  to  determine  the  mean  velocity 
relations  (Keulegan  1938).  For  wide  channels,  with  essentially  two- 
dimensional  flow,  the  mean  velocity  relation  is 

V  2.3  ,  11. ID  „  „ 

«-  '  —  108  "T—  2'7 
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where  V  is  the  average  flow  velocity. 

Several  difficulties  arise  in  application  of  the  logarithmic  veloc¬ 
ity  laws  to  rough  surfaces. 

1.  Origin  for  Velocity  Profile.  Several  investigators,  including 
Einstein  and  El-Samni  (1949),  O'Loughlin  and  McDonald  (1964), 
Knight  and  McDonald  (1979),  Bayazit  (1982),  and  Coleman,  Hodge, 
and  Taylor  (1984),  have  shown  that  the  velocity  profile  origin 
for  rough  surfaces  lies  below  the  tops  of  the  roughness  ele¬ 
ments.  There  is  no  general  agreement  as  to  the  location  of  the 
origin.  The  relation  between  velocity  and  tractive  force  is 


'.VJWvVV 
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sensitive  to  the  origin  location,  particularly  at  high  relative 


roughness . 


2.  K  Value.  Previous  studies  have  used  K  values  ranging  from 
s  s 


d^Q  (OCE  1970)  to  3.5dg^  (Hey  1979).  Particle  sizes  d^^  , 


dR,  ,  etc.,  refer  to  the  size  of  which  a  given  percentage  is 


finer  by  weight.  Kamphius  (1974)  found  Kg  =  2dg^  for  depth/ 


dQn  >  10  .  Van  Rijn  (1982)  determined  an  average  value  of 


Ks  =  3d90  ‘ 


3.  Effects  of  Relative  Roughness.  Yalin  (1977)  has  shown  that 


Equation  2.6  is  not  valid  at  relative  depth  D/d^ ^  less  than 


approximately  10  because  Ks/dgo  varies  below  D/d^  =  10  • 


Other  investigators  have  also  suggested  limiting  application  of 


the  logarithmic  velocity  equations  to  small  scale  roughness. 


Bathurst,  Graf,  and  Cao  (1982)  give  D/dg^  >  6  for  small-scale 


roughness.  Van  Rijn  (1982)  places  the  strictest  requirement  by 


limiting  application  of  the  logarithmic  velocity  laws  to 


D/Kg  >  10  .  Van  Rijn  (1982)  found  Kg  =  3d^^  which  implies  a 


limitation  D/dgQ  >  30  on  the  logarithmic  laws. 


Von  Karman  <  .  There  has  been  considerable  disagreement  over 


the  von  Karman  k  and  its  constancy  in  clear  versus  sediment¬ 


laden  flow.  Coleman  (1981)  found  that  by  evaluating  k  in  the 


lower  15  percent  of  the  flow,  k  was  the  commonly  used  0.4  for 


clear  or  sediment-laden  flow.  However  for  high  relative  rough¬ 


ness  D/d  =  4.0  and  8.5,  Bayazit  (1982)  found  <  <  0.4  for 


clear  water  flow  in  the  region  near  the  bed.  Uram  (1981)  found 


von  Karman’ s  <  both  higher  and  lower  than  0.4  depending  on 


the  nature  of  the  roughness. 


Summarizing,  other  investigators  have  suggested  that  a  constant 

Shields  coefficient  and  the  logarithmic  velocity  laws  should  not  be  used 

for  problems  involving  high  relative  roughness. 

2.2  EXISTING  CRITICAL  VELOCITY  METHODS 
FOR  PARTICLE  STABILITY 

Some  of  the  earliest  stability  relations  used  particle  size  or 
weight  as  a  function  of  velocity.  Graf  (1971)  presented  the  general 
relation 

2 

2Kg  (tan  $  cos  a  -  sin  a) 

Cd  Kj  +  Ct  K2  tan  T  2'8 

where 

=  bottom  velocity 
pg  =  stone  density 
=  coefficients 

a  =  bottom  angle  with  horizontal  in  flow  direction 
*  drag  coefficient 
C^  =  lift  coefficient 

Graf  referred  to  the  right  side  of  Equation  2.8  as  the  sediment  coeffi¬ 
cient  which  varies  with  particle  characteristics  (shape,  size,  unifor¬ 
mity,  texture,  repose  angle)  and  flow  characteristics. 

Forchheimer  (1914)  reported  that  as  early  as  1753,  A.  Brahms 
presented  the  relation 

V,  =  C_  W1/6  2.9 

D  J 

where  W  is  the  unsubmerged  stone  weight.  Equation  2.9  is  a  simple 
form  of  Equation  2.8.  Isbash  (1935)  related  stone  size  for  dam  closures 
to  a  bottom  velocity  called  the  "velocity  against  the  stone." 


Equation  2.8  is  the  form  used  by  Isbash  and  serves  as  the  basis  for 
Hydraulic  Design  Criteria  (HDC)  Sheet  712-1  (US  Army  Corps  of  Engineers). 
Average  velocity  is  used  in  HDC  712-1  instead  of  bottom  velocity,  which 
may  cause  these  curves  to  be  rather  conservative  for  low  turbulence 
flows.  The  National  Crushed  Stone  Association  (1978)  presents  guidance 
for  sizing  riprap  based  on  average  velocity.  The  California  Division  of 
Highways  (1970)  uses  a  design  equation  having  the  same  form  as 
Equation  2.8. 

Blodgett  and  McConaughy  (1986)  proposed  the  following  relation  for 
stable  rock  size  based  on  extensive  prototype  data 

2  44 

dc~  =  0.01  V  '  2.10 

50  a 

where  V  is  the  cross-section  average  channel  velocity.  Adjustments 
Si 

for  bank  angle,  unit  stone  weight,  channel  shape,  etc.,  are  not  used  in 
this  design  procedure. 

Critical  velocity  relations  using  average  velocity  and  depth  are 
also  used  for  particle  stability.  They  have  been  rewritten  in  a  common 
form  to  assist  in  their  comparison.  Straub  (1953)  presented  the  average 
velocity  and  depth  relation 


Neill  (1967)  used  dimensional  analysis  to  determine  the  pertinent  rela¬ 
tionships  for  stability  of  coarse,  uniform  bed  material  and  conducted 
scour  tests  using  the  incipient  criterion  of  first  movement  by  visual 
observation.  His  conservative  design  curve  is  represented  by  the 


equation 


(  y  w  ^ 

|1/2  V 

Vs  - 
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Neill  and  Van  Der  Giessen  (1966)  suggested  that  relative  roughness, 
which  results  from  the  dimensional  analysis,  is  connected  with  the 
intensity  of  turbulent  fluctuation.  Neill  (1968)  stated  that  because 
the  flume  size  and  test  section  area  were  constant,  the  first  movement 
criterion  was  more  severe  for  the  smaller  particles  and  Equation  2.12 
may  not  be  valid.  Because  the  test  section  contained  smaller  particles, 
and  therefore  more  particles,  a  greater  probability  of  movement  exists. 
Neill  (1968)  also  stated  that  the  equation  is  applicable  to  problems 
such  as  riprap  stability.  Bogardi  (1968)  presented  particle  stability 
data  covering  a  wide  range  of  d/D  and  determined  the  relation 


5  =  °-26 


which  is  almost  identical  to  Neill  (1967).  Cooper  (1970)  analyzed 
sediment  transport  data  for  low  rates  of  transport  (concentration  = 

1  part  per  million)  and  found  good  agreement  with  Neill's  (1967)  rela¬ 
tion.  Grace,  Calhoun,  and  Brown  (1973),  Maynord  (1978),  and  Reese 
(1984)  used  the  riprap  stability  relation 


Y  -  Y 
s  w 


which  is  identical  to  Straub  (1953). 

Combining  and  rearranging  Equations  2.4  and  2.7  results  in  the  OCE 
(1970)  procedure  for  riprap  design  using  average  velocity  and  depth: 
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With  the  appropriate  coefficients,  Equations  2.14  and  2.15  give  similar 
results  over  a  wide  range  of  d/D  .  Reese  (1984)  demonstrated  that 
these  two  relations  differ  only  by  the  velocity  profile  used.  Equa¬ 
tion  2.14  is  based  on  a  power  velocity  profile  while  Equation  2.15  is 
based  on  a  logarithmic  velocity  profile. 

Determining  which  velocity  to  use  is  an  important  step  in  develop¬ 
ing  a  riprap  sizing  method  based  on  velocity.  Some  form  of  bottom 
velocity  is  the  most  representative  because  it  is  closest  to  the  bed. 
However  bottom  velocities  are  difficult  to  predict  and  measure  (Bogardi 
1978)  because  the  velocity  near  the  bottom  varies  rapidly  with  distance 
from  the  bed.  Surface  velocities  are  easy  to  measure  but  difficult  to 
predict  and  are  not  representative  because  they  are  far  removed  from  the 
bed.  Bogardi  (1978)  recommended  the  use  of  mean  velocity  in  critical 
velocity  relations.  Mean  velocity  is  the  easiest  to  calculate  using 
both  numerical  and  physical  modeling  techniques. 


2.3  PREVIOUS  STUDIES  ON  THE  EFFECTS 
OF  GRADATION,  THICKNESS,  AND 
SHAPE  ON  RIPRAP  STABILITY 


The  effects  of  gradation  on  particle  stability  or  resistance  are 
generally  accounted  for  by  determining  a  characteristic  size  which  repre¬ 
sents  any  gradation.  In  the  case  of  resistance,  the  larger  size  frac¬ 
tions  are  generally  used  for  the  characteristic  size  (van  Rijn  1982, 
Bayazit  1982).  In  the  case  of  stability,  the  characteristic  size  is 


found  to  vary.  Einstein  (1942)  found  d^  to  be  the  effective  size  for 


movement  of  sand  mixtures.  Schoklitsch  (1962)  used  d^^  in  stability 


.VI 


relations.  The  California  Division  of  Highways  (1970)  used  in  the 

riprap  sizing  relation.  Peterka  (1958)  used  d^^  in  the  riprap  sizing 
relation  for  placement  downstream  of  stilling  basins.  Shen  and  Lu 
(1983)  found  d^Q  to  be  the  characteristic  size  of  nonuniform  surface 
material  on  an  armored  bed.  Shen  and  Lu  suggested  that  increaced  turbu¬ 
lence  caused  by  the  larger  particles  decreases  the  stability  of  non- 
uniform  materials.  Anderson,  Paintal,  and  Davenport  (1968)  conducted 
flume  tests  showing  that  nonuniform  ripraps  are  less  stable  than  uniform 
ripraps  having  the  same  average  size.  These  results  show  that  the  char¬ 
acteristic  size  is  less  than  the  average  size.  Maynord  compared  the 
stability  of  various  riprap  gradations  and  found  that  d^^  was  the 
characteristic  size  for  riprap  placed  to  a  thickness  of  ld^^  •  How¬ 
ever,  these  tests  differed  from  prototype  placement  of  riprap  because 
the  careful  placement  techniques  used  in  the  model  prevented  segregation 
of  sizes  with  the  nonuniform  ripraps.  Many  riprap  sizing  relations  have 
used  as  the  characteristic  size  (OCE  1970,  Anderson,  Paintal,  and 

Davenport  1968,  US  Department  of  Transportation  1975,  Blodgett  and 
McConaughy  1986) . 

Standardized  riprap  gradations  have  been  used  by  OCE  (1971), 
California  Division  of  Highways  (1970),  and  the  US  Army  Engineer  Divi¬ 
sion,  Lower  Mississippi  Valley  (1982).  Simons  and  Senturk  (1977)  and 
the  US  Department  of  Transportation  (1975)  present  a  single  curve 
defining  riprap  gradation. 

Studies  were  not  found  on  the  effects  of  varying  blanket  thickness 
on  riprap  stability.  Present  OCE  (1971)  guidance  requires  a  thickness 
of  IdjQQ  (maximum)  or  1.58,.^  (maximum),  whichever  is  larger,  for 


placement  in  the  dry. 
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Shape  effects  on  riprap  stability  are  important  in  determining 
which  shapes  are  acceptable.  Neill  (1968)  compared  the  stability  of 
spheres  and  "irregular  grains"  and  found  no  significant  difference  if 
the  equivalent  spherical  diameter  (volume  basis)  was  adopted  as  the  size 
of  the  irregular  grains.  Olivier  (1967)  conducted  tests  on  overflow 
rock  dams  and  found  that  rounded  stone  had  to  be  approximately  15  per¬ 
cent  larger  than  crushed  stone  for  equivalent  stability.  This  was 
attributed  to  surface  smoothness,  not  shape.  Present  OCE  (1970) 
guidance  for  riprap  shape  is  as  follows: 

1.  Stone  predominantly  angular 

2.  No  more  than  25  percent  of  stones  having  a  stone  length  l  to 
stone  thickness  b  ratio  of  >  2.5 

3.  No  stone  having  H/b  >  3.0 

2.4  EFFECTS  OF  SIDE  SLOPE  ON 

PARTICLE  STABILITY 

Since  most  riprap  is  placed  on  channel  banks,  the  influence  of  side 
slope  angle  on  riprap  stability  is  important.  Carter,  Carlson,  and  Lane 
(1953)  presented  the  effects  of  side  slopes  on  particle  stability  by 
defining  forces  parallel  and  normal  to  the  angle  of  repose  of  the  mate¬ 
rial.  The  equilibrium  condition  given  by  Carter,  Carlson,  and  Lane  is 


where 

Wg  =  submerged  weight  of  stone 
0  =  angle  of  side  slope  with  horizontal 
a  =  effective  area  of  particle 


tan  $  = 


,  2  2  2  2 
W  Z  sin  9  +  a  T  * 
/  J s _ s_ 

W  cos  0 


Tg  =  critical  tractive  force  for  particle  on  side  slope 
Carter,  Carlson,  and  Lane  defined  the  tractive  force  ratio  K  as  the 
ratio  of  force  on  sloping  side  to  that  on  level  surface  necessary  to 
cause  impending  motion 


K 


tan^  0 
2 

tan  $ 


sin^  Q 
sin^  $ 


2.17 


Equation  2.17  is  used  in  many  riprap  design  procedures  including 
Anderson,  Paintal,  and  Davenport  (1968),  US  Department  of  Transportation 
(1975),  and  OCE  (1970,  1971). 

An  alternate  formulation  by  Graf  (1971)  includes  lift  force 
and  the  angle  of  inclination  of  the  drag  force  or  shear  stress  as  a 
result  of  secondary  motion  8  ,  which  is  especially  pronounced  in 
channel  bends.  The  equilibrium  condition  is  alternately  written 


tan 


$  = 


sin^  0  +  2aT  W  sin  0  sin  B  +  a^T  ^ 
s  s  s  s 


W  cos  0  -  Ft 
s  L 


2.18 


Lack  of  information  on  the  angle  B  has  prevented  evaluation  of  this 
form  of  the  side  slope  stability  analysis.  Christensen  (1972)  developed 
a  side  slope  stability  analysis  which  included  lift  and  showed  that  the 
relation  given  by  Equation  2.17  is  not  conservative.  Stevens  and  Simons 
(1971)  determined  the  stability  of  coarse  particles  on  a  side  slope 
based  on  equilibrium  of  moments  instead  of  forces.  Relative  safety 
factors  can  be  determined  with  this  method  and  the  authors  concluded 


that  the  Carter,  Carlson,  and  Lane  (1953)  method  yields  larger  sizes 
than  required  by  the  Stevens  and  Simons  method. 


No  investigations  were  found  that  test  these  side  slope  equations 


for  open  channel  flow.  There  have  been  tests  in  the  wave  environment 
that  test  the  applicability  of  the  Hudson  (1958)  equation,  which 
follows: 


W  = 


2.19 


where 


H  =  wave  height 
Kp  =  stability  coefficient 

Since  wave  forces  act  up  and  down  the  side  slope,  the  effects  of  side 
slope  angle  are  expected  to  be  more  severe  than  that  in  open  channel 
flow  where  forces  act  along  the  slope.  Comparing  1V:1.5H  and  1V:3H  side 
slopes  in  Hudson's  equation  gives  the  wave  effect: 


d  (IV: 1 ,5H) 
d  (IV: 3H) 


1.26 


Using  the  open  channel  Equation  2.17  with  0  =  40  deg 
the  velocity  effect: 


(OCE  1970)  gives 


d  ( IV: 1 . 5H) 
d  (1V:3F) 


This  comparison  suggests  that  the  tractive  force  relation  (which  has  not 
been  tested  against  stability  data)  overestimates  the  effects  of  side 
slope  angle  on  stability. 

2 . 5  SUMMARY 


Several  investigators  have  proposed  that  Shields  coefficient  should 


A  table  of  factors  for  converting  non-SI  units  of  measurement  to 
SI  (metric)  units  is  found  on  page  iii. 


vary  with  relative  roughness.  Many  of  the  existing  riprap  design  tech¬ 
niques  use  a  constant  Shields  coefficient. 

Past  studies  have  shown  that  the  logarithmic  velocity  laws  should 
be  limited  to  small-scale  roughness.  Present  riprap  guidance  does  not 
place  any  limitations  on  use  of  these  laws.  Other  factors,  including 
determining  the  correct  values  of  Kg  ,  ic  ,  and  the  profile  origin, 
compound  the  difficulty  in  using  these  laws  for  surfaces  having  high 
relative  roughness. 

Several  different  velocity-based  riprap  sizing  methods  have  been 
developed.  Average  velocity  is  recommended  for  use  in  these  equations. 

Previous  studies  on  gradation  effects  on  the  stability  of  riprap 
have  used  a  characteristic  size  ranging  from  d^  to  d^  .  No  studies 
were  found  addressing  the  effects  of  riprap  thickness  on  stability. 

Side  slope  stability  equations  have  used  equilibrium  of  both  forces 
and  moments.  Information  was  not  found  in  which  these  equations  were 
tested  against  stability  data.  Comparison  of  the  side  slope  equations 
for  open  channel  flow  with  equations  tested  against  wave  data  suggests 
that  the  existing  side  slope  relations  for  open  channel  flow  over¬ 
estimate  the  effects  of  side  slope  angle. 


CHAPTER  3 


EXPERIMENTAL  INVESTIGATION 

Experimental  studies  were  conducted  to  determine  the  stability  and 
resistance  to  flow  of  riprap  having  gradation,  thickness,  and  shape 
similar  to  riprap  used  in  the  prototype  installations.  This  chapter 
describes  facilities,  model  riprap,  failure  criteria,  test  procedures 
and  data  collection,  and  data  restrictions.  Additional  information  on 
the  Colorado  State  University  (CSU)  studies  can  be  found  in  Fiuzat, 

Chen,  and  Simons  (1982),  Fiuzat  and  Richardson  (1983),  Ruff  et  al. 
(1985),  and  Ruff  et  al.  (1987). 

3.1  TEST  FACILITIES 

One  flume  at  CSU,  Fort  Collins,  Colorado,  and  three  flumes  at  the 
US  Army  Engineer  Waterways  Experiment  Station  (WES),  Vicksburg, 
Mississippi,  were  used  to  conduct  the  riprap  tests.  The  CSU  flume  is 
200  ft  long  by  8  ft  wide  by  4  ft  deep  and  can  be  tilted  from  0  to  2  per¬ 
cent  bottom  slope.  Maximum  discharge  is  100  cfs.  The  sides  and  bottom 
of  the  flume  are  made  primarily  of  aluminum.  A  portion  of  the  side  of 
the  flume  is  made  of  Plexiglas  to  allow  observation  of  the  test  section. 
Two  gates  installed  at  the  downstream  end  of  the  flume  allow  control  of 
the  water  level  in  the  flume  under  subcritical  flow  conditions.  A 
motorized  carriage  can  travel  along  the  flume  for  carrying  data  collec¬ 
tion  instruments  and  photographic  equipment.  A  schematic  diagram  of  the 
flume  and  the  test  section  is  shown  in  Figure  3.1.  The  initial  100  ft 


Flow  Development  Section  Transition  Test  Section 

Hesd 


of  the  flume  was  used  for  flow  development  and  transition  into  the 
test  section. 

The  CSU  tests  consisted  of  four  phases.  Phases  I-III  addressed 
stability  of  bottom  riprap  having  varying  gradation,  thickness,  and 
shape.  The  Phase  IV  tests  addressed  stability  of  side  slope  riprap.  In 
Phase  I  and  II  test  series,  large  6-  to  10-in.  rocks  cemented  to  the 
flume  floor  between  stations  0  and  80  produced  a  fully  developed  hydrau¬ 
lically  rough  boundary  flow  at  the  beginning  of  the  20-ft  transition. 
Rock  similar  in  size  to  that  in  the  test  section  was  placed  in  the  20-ft 
transition  to  eliminate  the  abrupt  change  In  roughness  between  the  flow 
development  section  and  the  test  section.  In  the  Phase  III  test  series, 
the  large  6-  to  10-in.  rocks  were  placed  in  the  initial  60-70  ft  of  the 
flume.  A  40-ft-long  transition  was  used  in  the  Phase  III  tests.  The 
test  section  varied  from  40  to  50  ft  in  length  for  Phases  I-III.  Detail 
of  the  Phase  IV  test  facility,  in  which  a  IV: 2H  side  slope  was  tested, 
are  shown  in  Figure  3.2. 

The  WES  trapezoidal  channel  model  is  described  in  Maynord  (1978). 


This  facility  had  a  5-ft  bottom  width  with  1V:2H,  1V:3H,  and  1V:4H  side 
slopes.  Discharge  capacity  was  35  cfs,  and  a  constant  bottom  slope  of 
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Figure  3.2.  CSU  Phase  IV  side  slope  test  flume 
0.008  ft/ft  was  used  in  all  tests.  A  tailgate  was  used  to  control  depth 
of  flow. 

The  WES  tilting  flume  is  3  ft  wide  by  1  ft  deep  by  75  ft  long. 
Maximum  discharge  is  5.6  cfs.  Bottom  slope  can  be  varied  from  0  to 
2.2  percent,  and  a  tailgate  at  the  downstream  end  of  the  model  is  used 
to  control  depth  of  flow  for  subcritical  flows.  Steel  rails  set  to 
grade  are  used  to  support  instrumentation  devices. 

The  WES  curved  channel  model  is  shown  in  Figure  3.3.  This  trape¬ 
zoidal  channel  has  two  100-deg  bends  with  a  centerline  radius  of  22  ft. 
The  bends  are  separated  by  a  15-ft  straight  reach,  and  the  straight 
reach  on  each  end  of  the  channel  is  25  ft  in  length.  The  bottom  width 
is  7.0  ft,  and  side  slopes  are  1V:2H.  The  bottom  slope  is  0.0025  ft/ft, 
and  the  discharge  varies  up  to  15  cfs. 
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Figure  3.3.  WES  curved  channel  model 

3.2  MODEL  RIPRAP 

The  characteristics  of  the  model  riprap  used  in  these  investiga¬ 
tions  are  given  in  Table  3.1.  Gradations  for  the  CSU  flume  are  shown  in 
Figures  3. 4-3. 7. 

All  model  riprap  was  crushed  rock.  Shape  characteristics  of  the  model 
riprap  are  shown  in  Table  3.2. 

3.3  FAILURE  CRITERIA 

At  the  outset  of  these  experimental  studies,  an  acceptable  failure 
criterion  had  to  be  determined.  The  selected  failure  criterion  must  be 
able  to  be  used  to  determine  riprap  stability  for  a  range  of  riprap 
gradation  and  blanket  thicknesses.  Most  sediment  transport  studies  using 
uniform  materials  have  weighed  the  transported  material  for  various  flow 
rates  and  extrapolated  the  transport  rate  to  zero  to  determine  what  is 


Figure  3.6.  Size  distribution  for  CSU  Phase  III 
termed  "incipient  motion."  Applying  this  technique  to  different  riprap 
blanket  thicknesses  would  probably  yield  little  variation  with  thick¬ 
ness.  Applying  this  technique  to  nonuniform  ripraps  would  give  biased 
results  because  some  of  the  finer  material  in  nonuniform  ripraps  will  be 
moved  without  ultimate  failure  of  the  riprap  revetment. 

Another  existing  failure  criterion  is  the  technique  used  by  Neill 
(1967),  which  was  a  visual  observation  of  first  movement.  This  tech¬ 
nique  would  be  successful  for  uniform  materials  but  unsuccessful  for 
nonuniform  (graded)  ripraps  of  varying  thickness.  The  idea  of  painting 
rocks  in  the  test  section  was  rejected  because  it  would  yield  no  infor¬ 
mation  about  the  effects  of  thickness  for  nonuniform  ripraps. 


Figure  3.7.  Size  distribution  for  CSU  Phase  IV 
An  important  consideration  in  riprap  stability  is  that  the  under¬ 
lying  material  should  not  be  exposed  to  the  forces  of  the  flowing  water. 
It  is  not  important  if  some  of  the  finer  material  resting  on  the  surface 
of  a  nonuniform  riprap  is  washed  away.  Another  factor  which  must  be 
considered  with  riprap  stability  is  size  segregation  during  placement. 
The  selected  failure  criterion  must  be  able  to  address  the  effects  of 
size  segregation  when  using  nonuniform  materials. 

To  meet  these  requirements,  the  concept  of  incipient  failure  is 


used  in  this  investigation  to  define  the  flow  conditions  at  which  any 
portion  of  the  underlying  material  has  been  exposed.  Use  of  this 
failure  criterion  allows  determination  of  the  stability  of  various 


gradations  and  thicknesses.  It  is  the  only  failure  criterion  which  was 
considered  to  address  the  effects  of  size  segregation.  The  incipient 
failure  criterion  is  not  the  same  as  the  incipient  motion  criterion  used 
in  sediment  transport  studies.  Incipient  failure  defined  the  flow  con¬ 
ditions  which  lead  to  failure  of  the  riprap  blanket.  Incipient  motion 
defines  the  flow  conditions  at  which  the  rate  of  particle  movement 
approaches  zero.  Incipient  motion  could  not  be  used  in  this  study 
because  it  would  not  allow  determination  of  the  effects  of  riprap  grada¬ 
tion  or  thickness. 

3.4  TEST  PROCEDURE  AND  DATA  COLLECTION 

In  the  CSU  and  WES  trapezoidal  channel  tests,  a  fabric  was  used  to 
separate  the  riprap  from  the  bed  of  underlying  sand.  In  the  WES  tilting 
flume  tests,  the  fabric  was  placed  directly  on  the  flume  floor.  While 
the  riprap  was  being  placed,  the  riprap  surface  was  not  tamped  or  packed 
to  best  simulate  prototype  placement.  The  flow  conditions  at  which  the 
rock  would  fail  were  estimated  using  existing  riprap  sizing  techniques. 
The  initial  test  began  with  low  flow  rates  and  slopes  well  below  the 
estimated  failure  condition.  The  riprap  was  tested  for  2  hours,  after 
which  the  test  section  was  examined  for  any  exposed  areas  of  the  under¬ 
lying  fabric.  If  no  exposure  of  the  fabric  occurred,  the  flow  rates  or 
slope  was  increased  and  the  2-hour  test  repeated.  This  process  was 
repeated  until  the  fabric  was  exposed.  After  the  test  section  was 
repaired,  the  previous  stable  slope  was  run  for  4-8  hours  to  ensure  sta¬ 
bility  of  the  riprap.  In  case  of  failure,  the  slope  and/or  discharge 
was  further  reduced  and  another  4-  to  8-hour  run  was  conducted  until 
stable  conditions  were  found.  The  WES  tilting  flume  tests  differed  in 


that  near  the  point  of  failure,  all  runs  used  iu  the  analysis  were  4-  to 
8-hour  runs. 

In  the  CSU  and  WES  tilting  flume  studies,  uniform  flow  was  main¬ 
tained  by  adjusting  the  tailgate  at  the  downstream  end  of  the  flume  for 
subcritical  flows.  The  WES  trapezoidal  channel  tests  had  a  mild,  grad¬ 
ually  varied  flow  regime  because  of  the  lack  of  slope  variability.  Flow 
uniformity  in  the  WES  curved  channel  model  was  maintained  by  keeping  the 
same  depth  at  the  upstream  and  downstream  ends  of  the  model. 

During  the  tests  at  both  CSU  and  WES,  discharges  were  measured  by 
calibrated  venturi  and  orifice  meters,  velocity  was  determined  with 
pitot  tubes  and  propeller  meters,  and  depths  were  measured  with  point 
gages. 

In  the  CSU  tests,  a  "general  datum"  for  each  rock  thickness  was 
established  by  the  following  procedures: 

1.  The  flume  was  set  to  the  horizontal  position. 

2.  Water  was  added  to  the  flume  until  about  90  percent  of  the 
rocks  were  covered  with  water. 

3.  The  elevation  of  the  water  surface  was  measured  at  the  loca¬ 
tions  where  flow  depths  were  measured. 

4.  These  elevations  were  considered  as  the  elevations  of  the  bot¬ 
tom  of  channel  (general  datum)  in  measuring  the  flow  depth. 

In  the  WES  tests,  the  datum  was  set  by  placing  a  flat  plate  of  known 
thickness  on  top  of  the  riprap  surface  to  establish  the  datum. 

3 . 5  DATA  RESTRICTIONS 

Two  areas  of  concern  generally  surface  in  the  course  of  any  flume 
investigation.  First,  flow  conditions  must  be  turbulent  to  ensure  that 
viscous  forces  are  insignificant  in  the  flume  just  as  they  are  in  the 
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prototype.  To  ensure  rough  turbulent  flow,  the  following  restrictions 
were  placed  on  the  data  to  be  used  in  the  analysis: 
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>  2.5(10) 


(O'Loughlin  et  al.  1970) 


The  second  area  of  concern  in  flume  studies  is  how  to  handle  the 
effects  of  the  sidewalls.  Previous  sediment  transport  studies  have 
frequently  used  the  sidewall  correction  procedure  given  in  Vanoni 
(1977).  This  sidewall  correction  procedure  results  in  the  average  bed 
shear  stress.  This  method  takes  the  central  region  of  the  flume  (where 
the  flow  is  essentially  two-dimensional)  and  the  two  side  regions  (where 
the  shear  stress  and  velocities  are  reduced)  and  determines  a  weighted 
average.  In  this  type  of  study,  the  riprap  generally  fails  and  veloci¬ 
ties  are  measured  in  the  central  region  of  the  flume.  What  is  needed  is 
not  the  weighted  average  but  the  values  of  shear  stress  and  velocity  in 
the  central  portion  of  the  flume.  The  velocities  pose  no  difficulty 
because  they  are  measured  in  the  central  portion  of  the  flume,  but  shear 
stress  needs  to  be  calculated.  If  the  central  portion  of  the  flume  is 


sufficiently  wide,  then  the  shear  stress  is  best  approximated  by  y  DS  . 

w 


To  ensure  that  the  central  region  is  wide  enough,  Neill  (1967)  and  van 
Rijn  (1982)  required  that  the  ratio  of  flume  bottom  width  B  to  depth 
(aspect  ratio)  be  equal  to  or  greater  than  5.  As  part  of  this  study, 
the  limiting  aspect  ratio  was  evaluated  with  velocity  measurements  taken 
in  a  straight,  riprapped  bottom,  smooth  sidewall  flume.  Detailed 
velocity  measurements  were  taken  at  aspect  ratios  of  4.0,  4.9,  and  7.3 
(Figure  3.8).  The  tests  with  aspect  ratios  of  4.9  and  7.3  show  a 
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1.  The  CSU  tests  were  generally  conducted  with  sequential  dis¬ 
charges  of  25,  50,  75,  and  100  cfs.  By  the  time  the  75-  and 
100-cfs  tests  were  conducted,  the  riprap  in  the  flume  was 
"well-seasoned."  Any  weak  spots  had  already  failed  and  had 
been  repaired.  The  tests  conducted  at  these  higher  discharges 
generally  did  not  meet  the  width/depth  >5  restriction. 

2.  At  the  deeper  depths,  slopes  were  mild  at  the  point  of  failure 
of  the  riprap.  Only  three  depths  were  measured  along  the  test 
section  for  each  test,  which  made  it  difficult  either  to  assume 
that  the  bottom  slope  equaled  the  energy  slope  or  to  compute 
the  energy  slope.  At  mild  slopes,  errors  in  determining  energy 
slope  can  be  large.  At  steeper  slopes,  the  bottom  slope  domi¬ 
nates  the  energy  slope  and  errors  due  to  a  limited  number  of 
depth  measurements  are  small.  This  factor  was  probably  signif¬ 
icant  only  for  the  smaller  ripraps. 

Fortunately,  data  meeting  the  width/depth  >5  requirement  are  sufficient 
to  define  riprap  stability  for  the  majority  of  problems. 
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CHAPTER  4 

ANALYSIS  AND  RESULTS 

4.1  APPLICABILITY  OF  EXISTING  SIZING 
RELATIONS  USING  A  CONSTANT 
SHIELDS  COEFFICIENT  OR 
LOGARITHMIC  VELOCITY  LAWS 


The  review  of  previous  work  presented  in  Chapter  2  indicates  that 
numerous  investigations  have  proposed  that  Shields  coefficient  should 
vary  with  relative  roughness.  Experimental  results  from  the  WES  and  CSU 
tilting  flumes  were  used  to  evaluate  Shields  coefficient  as  a  function 
of  relative  roughness.  Results  from  the  WES  trapezoidal  channel  were 
not  used  because  the  test  section  was  not  long  enough  to  accurately  mea¬ 
sure  the  water-surface  slope  so  that  shear  stress  could  be  computed. 

Only  those  data  sets  covering  a  large  range  of  D/d  and  having  the  same 
thickness,  gradation,  and  shape  were  used  in  this  analysis.  Shields 
coefficients  computed  for  the  four  data  sets  meeting  these  requirements 
are  shown  in  Figures  4. 1-4.4.  The  data  used  in  Figures  4. 1-4.4  are 
listed  in  Tables  4. 1-4. 7.  Shields  coefficient  is  computed  using  a 
combination  of  Equations  2.1  and  2.4  or 
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V  50 

and  only  data  meeting  the  limitations  in  Chapter  3  are  used  in  the 
analysis.  In  a  comparison  of  these  results  to  the  Shields  (1936)  work, 
the  difference  in  stability  criteria  must  be  considered.  The  Shields 


(1936)  investigation  measured  low  rates  of  transport  and  extrapolated 
these  values  to  a  zero  rate  of  transport  to  obtain  incipient  conditions. 
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Figure  4.3.  Shields'  coefficient  versus  D/d50  , 
thickness  =  1.4d^Q  ,  ^85^15  =  2.1  and  2.3 

(data  from  Tables  4.4  and  4.5) 
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Figure  4.4.  Shields'  coefficient  versus  B/d^  > 
thickness  =  2.1d^Q  ,  dg^/d^  =  2.1  and  2.3 
(data  from  Tables  4.6  and  4.7) 
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This  investigation  used  the  incipient  failure  stability  criterion  given 
in  paragraph  3.3.  The  best-fit  lines  on  Figures  4. 1-4.4  are  drawn  to 
separate  failure  runs  from  stable  runs  and  are  not  the  result  of  regres¬ 
sion  techniques.  Three  of  the  four  data  sets  (Figures  4.1,  4.3,  and 
4.4)  show  a  significant  increase  in  Shields  coefficient  with  a  decrease 
in  D/d5Q  .  This  is  the  same  variation  proposed  by  several  investi¬ 
gators  cited  in  Chapter  2.  Over  the  range  of  D/d^  tested,  there  was 
no  indication  that  Shields  coefficient  approached  a  constant  value  as 
proposed  by  Bathurst,  Graf,  and  Cao  (1982)  and  Bettess  (1984).  The 
average  of  the  best-fit  lines  shown  in  Figures  4.1,  4.3,  and  4.4  show 
that  Shields  coefficient  should  vary  according  to 


which  can  be  compared  to  Meyer-Peter  and  Muller  (1948)  results  given  in 
Equation  2.5. 

Logarithmic  velocity  relations  are  used  in  riprap  design  to  relate 
velocity  to  shear  stress.  Several  references  cited  in  Chapter  2  show 
that  the  logarithmic  velocity  relations  are  not  applicable  to  high  rela¬ 
tive  roughness  and  should  be  limited  to  small-scale  roughness.  The  mean 
velocity  logarithmic  Equation  2.7  is  the  equation  most  frequently  used 
in  riprap  design  problems  and  will  be  evaluated  in  this  analysis.  The 
mean  velocity  relation  results  from  integration  of  the  point  velocity 
relation  over  the  entire  depth  of  flow.  This  is  one  problem  with  the 
mean  velocity  equations,  if  Coleman  (1981)  is  correct  in  saying  that  the 
point  velocity  logarithmic  equation  is  applicable  in  only  the  lower 
15  percent  of  the  depth.  Another  problem  is  that  the  origin  for  the 


velocity  profile  is  assumed  equal  to  the  tops  of  the  roughness  elements 


in  the  integration.  This  assumption  is  satisfactory  for  low  relative 
roughness  but  not  for  high  relative  roughness.  The  effects  of  both  of 


these  assumptions  are  lumped  into  the  determination  of  Kg  .  The 
experimental  data  collected  in  the  WES  and  CSU  flumes  were  used  to 
define  the  applicability  of  the  logarithmic  velocity  relations.  Analy¬ 
sis  of  Equation  2.7  was  similar  to  Yalin  (1977)  in  which  K^/d^^  is 
determined  as  a  function  of  relative  roughness.  Results  are  presented 
in  Figure  4.5  for  tests  with  no  movement  and  meeting  the  data  require¬ 
ments  given  in  Chapter  3.  Data  used  in  Figure  4.5  are  given  in 

Tables  4.2-4.12.  Results  show  that  K.  / d_ _  is  not  constant  over  the 

s  90 

range  of  data  used  in  this  investigation.  This  result  is  consistent 

with  Yalin' s  results  showing  the  point  velocity  logarithmic  relation 

(Equation  2.6)  inapplicable  for  D/d^  <  *0  . 

4.2  DEVELOPMENT  OF  CRITICAL 
VELOCITY  RELATION 

One  of  the  objectives  of  this  study  is  to  develop  a  riprap  sizing 
method  based  on  velocity.  Dimensional  analysis  is  used  to  define  the 
dimensionless  variables  based  on  the  selection  of  all  relevant  param¬ 
eters.  The  dimensional  analysis  is  similar  to  that  proposed  by  Neill 
(1967)  in  which  mean  velocity  is  used  instead  of  the  critical  tractive 
force  approach  used  by  many  investigators.  The  relevant  parameters 
governing  the  stability  of  riprap  in  open  channels  are  as  follows: 
d  =  characteristic  particle  size,  L 
D  =  flow  depth,  L 
p  =  fluid  density,  M/L^ 

3 

Pg  =  stone  density,  M/L 

V  =  mean  velocity,  L/T 

V  =  absolute  viscosity,  M/LT 
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Figure  4.5.  Variation  of  K  / d--  with  relative  depth 

s  9U 

(data  from  Tables  4.2-4.12) 
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g  =  gravitational  acceleration,  L/T 
=  side  slope  factor 

SIDE 

S  =  channel  slope 
dQ  /d._  =  gradation  uniformity 

OJ  1 J 

N  =  blanket  thickness/d^Q 
F  =  shape  factor  and  surface  texture 

bHAr  h 

M,L,T  =  fundamental  dimensions  of  mass,  length,  and  time, 
respectively 

The  mean  velocity  in  this  investigation  is  the  average  velocity  in  the 
vertical  at  the  point  of  interest.  With  this  concept,  the  effects  of 
channel  alignment  or  curvature  can  be  incorporated  into  the  design  pro¬ 
cedure.  The  designer  must  determine  the  velocity  at  the  point  of  inter¬ 
est,  not  average  cross-section  values,  in  order  to  determine  rock  size. 
Methods  to  determine  this  average  velocity  at  the  point  of  interest  are 
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presented  in  Chapter  5.  Neither  g  nor  pg  can  be  independent  param¬ 
eters;  they  must  occur  in  the  combination  g(p  -  p)  which  is  the  sub- 

s 

merged  specific  weight  of  the  riprap  .  Replacing  p  with  v  = 
p/p  ,  the  relevant  parameters  can  be  written 


/(d,  D,  p,  V,  Yg »  v,  Fside.  S,  dg5/d15>  N,  FSHAPE)  =  0  4.3 

Out  of  these  eleven  variables  there  are  six  dimensional  variables  (d, 

D,  P,  V,  y  ' «  v)  and  five  dimensionless  quantities  (F.,-.^.,,  S,  doc/d  c,  N, 
s  blue  oj  ID 

F  ).  Since  there  are  three  fundamental  dimensions  (M,  L,  and  T) , 
b  HAr  c 

there  are  three  nondimensional  groups.  The  statement  can  be  rewritten 


/(V  *2’  *3*  FSIDE’  s*  d85^d15*  N’  FSHAPE^  °  4,4 


Using  repeating  variables  V  ,  D  ,  and  ,  the  tt  terms  are 
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Set  each  of  it's  equal  to  M°L°T°  and  solve  simultaneously  for  a  , 
b  ,  c  ,  and  d  .  This  results  in 
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'  The  statement  can  then  be  rewritten 
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!  The  Reynold's  number  term  VD/^  is  indicative  of  viscous  effects  which 

are  not  important  in  prototypes  and  in  the  model  sizes  used  in  this 
investigation.  The  influence  of  slope  is  important  for  steep  flows,  and 
Bathurst,  Graf,  and  Cao  (1982)  found  slopes  greater  than  2  percent  to 
have  significant  effects  on  incipient  conditions  of  bed  movement.  At 
the  condition  of  incipient  failure  of  riprap,  slope  and  particle 
size/depth  ratio  d/D  are  dependent.  A  steep  slope  implies  large  d/D 
at  incipient  conditions.  Since  this  investigation  is  limited  to  slopes 
equal  to  or  less  than  2  percent  and  since  d/D  is  retained  in  the 
analysis,  slope  is  omitted.  The  majority  of  open  channel  riprap 
problems  have  slopes  well  below  2  percent.  The  statement  of  relevant 
dimensionless  parameters  becomes 
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Riprap  stability  data  will  be  used  to  evaluate  the  importance  of  each  of 
these  parameters.  Channel  bottom  test  series  having  a  relatively  large 
range  of  d/D  and  having  the  same  gradation  uniformity  doc/d1c  ,  thick- 
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ness  N  ,  and  shape  F  were  used  to  evaluate 

bnAr  L 


d 

D 


gD 


function  of 


4.13 


Results  in  Figures  4. 6-4. 9  show  that  the  basic  equation  for  threshold  of 
incipient  failure  of  bottom  riprap  in  straight  channels  has  the  form 
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These  best-fit  lines  were  drawn  to  separate  stable  runs  from  failure 
runs  and  were  not  the  result  of  regression  techniques.  Equation  4.14  is 
the  same  form  found  by  Neill  (1967)  and  Bogardi  (1968)  and  will  be  used 
in  the  evaluation  of  the  effects  of  gradation,  thickness,  and  shape. 
Particle  size  d^  was  used  in  this  analysis  until  additional  analysis 
can  define  a  characteristic  size. 

An  equation  similar  to  Equation  4.14  can  be  derived  by  combining 
the  following  shear  or  tractive  force  relations: 


Manning's  equation 
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where  n  =  Manning's  resistance  coefficient  and  Strictler's  equation  is 
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Figure  4.6.  d^/D  versus  modified  Froude  number,  thickness 
doc/d. c  “  1.35  ,  WES  Phase  I  (Data  from  Table  4.1) 
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Figure  4.7.  d^/D  versus  modified  Froude  number,  thickness  =  Id^gg  , 
d  /d. c  =  2.8  and  2.5,  CSU  Phases  I  and  II  (Data  from  Tables  4.2  and 
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When  combined,  these  equations  yield 
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which  is  similar  to  Equation  4.14. 

Most  existing  riprap  design  procedures  fall  into  two  categories: 

1.  Constant  Shields  coefficient.  Equation  2.14  is  an  average 

velocity  relation  which  can  be  derived  using  a  constant  Shields 
coefficient  and  combining  equations  2.1,  2.4,  4.15,  4.16,  and 
4.17. 
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2.  Isbash  type  relations.  These  can  be  expressed  as 
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Many  US  Army  Corps  of  Engineers  (CE)  offices  have  charts 
relating  riprap  size  to  velocity  which  use  this  relation.  This 
relation  can  be  rewritten  in  the  form 
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using  the  full  form  of  the  Isbash  equation  and  dividing  both 
sides  by  depth. 

Comparing  Equations  2.14  and  4.20  to  the  equation  proposed  in  this 
investigation  (previously  proposed  by  Neill  (1967)  and  Bogardi  (1968)) 
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shows  that  riprap  stability  is  best  described  by  a  relationship  with  an 
exponent  that  falls  halfway  between  the  two  most  commonly  used  methods 
of  design. 

4.3  EFFECTS  OF  GRADATION,  THICKNESS , 

AND  SHAPE  ON  RIPRAP  STABILITY 

4.3.1  Gradation.  Variation  of  gradation  uniform-* t>  was  accounted 
for  by  using  a  characteristic  size  less  than  the  average  size  given  in 
several  references  cited  in  the  literature  review.  Size  segregation  can 
be  a  significant  factor  when  us-*ng  highly  nonuniform  materials  and  is 
probably  one  reason  the  characteristic  size  was  found  to  be  less  than 
the  average  size.  The  ratio  dQC/d,c  is  used  to  describe  the  unifor- 
mity  of  riprap  gradations.  Standard  CE  gradations  given  in  OCE  (1971) 
have  dg^/d^  =  1. 8-2.1  .  In  addition  to  the  results  presented  in 
Figures  4. 6-4. 9,  data  from  the  following  test  series  were  evaluated 


using  Equation  4.14 

(these  data 

sets  were  not 

used  in 

the  development 

Equation  4.14  because  they  do  not  cover  a  wide 

enough 

range  of  d/D) 

Source 

d85/d15 

Thickness 

Table 

Figure 

CSU  Phase  I 

3.9 

ldioo 

4.8 

4.10 

CSU  Phase  II 

4.6 

ld100 

4.9 

4.11 

To  evaluate  the  effects  of  gradation  for  riprap  placed  to  a  thick¬ 
ness  of  Icl^go  ’  coef f icionts  from  the  equations  shown  in  Fig¬ 

ures  4.6,  4.7,  4.10,  and  4.11  are  plotted  against  dg^/d^  in 
Figure  4.12.  Results  show  that  the  coefficient  varies  with  dg^/d^  , 
which  means  that  d^  is  not  the  characteristic  size  for  the  range  of 
gradations  tested.  Equation  4.14  was  evaluated  using  different 
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characteristic  sizes,  and  only  d 


(Figure  4.12)  was  shown  to  give  a 
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Figure  4.12.  Coefficient  in  Equation  4.14  versus 
^85/^15  *  Sickness  =  ld^^  ,  bottom  riprap 
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Figure  4.13.  d^/D  versus  modified  Froude  number, 
thickness  =  ld^  ,  d^/d^  =  1.35-4.6  ,  bottom 
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relatively  constant  C  in  Equation  4.14.  All  data  are  plotted  in 
Figure  4.13.  The  equation 
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4.21 


is  applicable  to  threshold  of  incipient  failure  of  riprap  placed  to  a 


thickness  of  ld^^  *  ^85^15  —  ’  ^q/D  =  0*020-0.25  ,  F  <  1.2  , 


on  the  bottom  of  straight  channels.  This  analysis  shows  that  either 


d^Q  can  be  used  in  Equation  4.14  with  a  coefficient  which  varies  with 
gradation  or  d^Q  can  be  used  in  Equation  4.21. 

4.3.2  Thickness.  Several  of  the  test  series  were  used  to  deter¬ 
mine  the  effects  of  blanket  thickness  on  riprap  stability.  Any  compari¬ 
son  of  different  thicknesses  of  riprap  must  be  conducted  with  the  same 
gradation.  Data  from  Maynord  (1978)  are  shown  in  Figure  4.14  for  a 


thickness  of  1.5d^Q  .  The  following  tabulation  summarizes  the  test 
series  used  in  the  analysis  of  thickness  effects: 


Source 

d85/d15 

Thickness 

d30/d50 

Table 

Figure 

WES  trapezoidal 
channel 

2.0 

1,5dioo 

0.83 

See  Maynord 
(1978) 

4.14 

CSU  Phase  III 

2. 1-2.3 

1,Adioo 

0.80 

4.4,  4.5 

4.8 

CSU  Phase  III 

2. 1-2.3 

2-ldioo 

0.80 

4.6,  4.7 

4.9 

The  coefficients  from  these  equations  are  determined  for  a  characteris¬ 
tic  size  d^Q  and  plotted  against  thickness  in  Figure  4.15.  For  thick¬ 
ness  of  l.OdjQQ  ,  the  coefficient  from  Equation  4,21  is  used.  Results 
show  that  increased  thickness  decreases  the  size  required  to  remain 
stable  up  to  a  thickness  of  2.0-2.5d^QQ  .  Additional  tests  are  needed 
to  evaluate  the  effects  of  thickness  for  other  gradations.  Note  that 
d^  was  shown  to  be  the  characteristic  size  for  a  thickness  of  l.Od^g 
only.  As  riprap  thickness  increases,  the  likelihood  of  areas  having 


Figure  4.14.  d^/D  versus  modified  Froude  number, 
thickness  =  1.5d^QQ  ,  dg^/d^  =  2.0  ,  WES  trapezoidal 


N 

Figure  4.15.  Coefficient  in  Equation  4.21  (using  d^)  ver 
blanket  thickness,  doc/d.c  =  2. 1-2.3  ,  bottom  riprap 


only  small  particles  (due  to  size  segregation)  decreases.  Another 
mechanism,  armoring,  may  also  exert  a  significant  influence  on  riprap 
stability.  These  may  be  the  reasons  that  thickness  is  seen  to  be  so 
significant  in  Figure  4.15. 

4.3.3  Shape .  Two  of  the  test  series  conducted  during  the  CSU 
Phase  III  tests  allow  a  comparison  of  the  effects  of  riprap  shape.  OCE 
(1970)  shape  guidance  requires  the  following: 

1.  Stone  predominantly  angular 

2.  No  more  than  25  percent  of  stones  having  l/b  >  2.5 

3.  No  stone  having  9jb  >  3.0 

Riprap  meeting  this  guidance  was  tested  and  compared  with  riprap  having 
the  following  characteristics: 

1.  All  stone  angular 

2.  Thirty  percent  of  stone  had  £/b  >  2.5 

3.  Eighteen  percent  of  stone  had  £ /b  >  3.0 

Results  of  these  two  test  series  are  plotted  in  Figure  4.16.  Data  used 

in  Figure  4.16  are  from  Tables  4.4  and  4.10.  Results  show  that  shape 

effects  are  insignificant  within  the  range  tested  in  this  investigation. 

Neill  (1967)  also  found  shape  effects  to  be  small.  The  stability  of 

rounded  rock  such  as  cobbles  was  not  addressed  in  this  investigation. 

4.4  EFFECTS  OF  SIDE  SLOPE  ON 
RIPRAP  STABILITY 

Three  areas  must  be  addressed  in  defining  the  effects  of  channel 
side  slopes  on  riprap  stability.  First,  the  effects  of  the  gravity  com¬ 
ponent  acting  downslope  and  the  influence  of  angle  of  repose  roust  be 
evaluated.  Second,  the  effects  of  the  side  slope  on  the  velocity  pro¬ 
file  and  distribution  must  be  incorporated  into  the  average  velocity 
relations  for  sizing  riprap.  Third,  side  slope  stability  tests  must  be 
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conducted  to  determine  the  combined  effect  of  the  gravity  component  and 
the  velocity  profile. 

4.4.1.  Effect  of  Gravity  Component  Acting  Downslope.  As  indicated 
in  the  literature  review,  several  different  methods  including  equilib¬ 
rium  of  moments  and  equilibrium  of  forces  have  been  used  to  define  the 
stability  of  a  particle  resting  on  a  channel  side  slope.  As  part  of 
this  investigation,  tests  were  conducted  in  the  WES  tilting  flume 
(Phase  IV)  to  compare  the  stability  of  riprap  resting  on  various  side 
slopes.  A  schematic  of  the  test  facility  is  shown  in  Figure  4.17.  The 
side  slope  was  hinged  at  the  bottom  of  the  slope  to  facilitate  changing 
the  side  slope.  The  riprap  surrounding  the  test  section  was  the  same 
size  as  used  in  the  test  section  and  was  glued  to  the  side  slope  to 
ensure  that  the  velocity  profile  and  turbulence  characteristics  of  the 
approach  flow  did  not  vary  from  test  to  test.  Results  for  six  different 
side  slopes  using  uniform  riprap  with  a  thickness  of  ld^^  are  s^own 
in  Table  4.13.  Bottom  velocity  was  used  to  define  the  imposed  velocity 
and  was  measured  2.9d^^  above  the  side  slope  as  shown  in  Figure  4.18. 
Results  show  decreasing  bottom  velocity  for  increasing  side  slope. 

The  tractive  force  ratio  as  used  by  Carter,  Carlson,  and  Lane 
(1953)  is 

K.Is.cose  /  i  -  HaL? ,  I,  .  (2.17  b±s) 

c  V  tan  $  \  sin  $ 

Given  the  same  fluid,  particle  characteristics,  and  depth,  shear  stress 
is  proportional  to  the  second  power  of  the  velocity 

T-CUV2  4.22 

The  WES  tilting  flume  side  slope  tests  were  conducted  with  the  same 
fluid,  particle  characteristics,  and  depth.  The  only  factor  that  varied 
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Figure  4.17.  WES  side  slope  test  flume 


Figure  4.18.  Location  of  bottom  velocity  measurements  in  WES 
tilting  flume  side  slope  tests 


was  the  side  slope  angle.  Having  established  these  conditions  Equa¬ 
tion  4.22  can  be  substituted  into  2.17  to  obtain 


K 


4.23 


where 

Vg  =  critical  velocity  for  particle  on  side  slopes 

Vc  =  critical  velocity  for  particle  on  horizontal  bed 
The  flattest  side  slope,  1V:4H,  and  the  horizontal  test  yield  essen¬ 
tially  the  same  critical  velocity  and  will  be  used  for  in  this 

analysis.  The  tractive  force  ratio  K  from  Equation  4.23  is  plotted 
against  the  side  slope  angle  0  in  Figure  4.19.  Also  shown  in  this 
figure  is  the  analysis  of  Carter,  Carlson,  and  Lane  (1953)  using  an 
angle  of  repose  of  40  deg  (OCE  1970).  The  Carter,  Carlson,  and  Lane 
method  shows  a  greater  decrease  in  stability  than  the  experimental  data. 
The  experimental  results  are  consistent  with  the  findings  of  Hughes, 
Urbonas,  and  Stevens  (1983)  stating  that  "rock  size  does  not  need  to  be 
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the  slope.  The  average  value  of  repose  angle  obtained  for  this  revet¬ 
ment  configuration  under  dry  conditions  was  52  deg  (see  Table  4.14). 
These  tests  were  repeated  with  the  test  section  submerged,  and  the  aver¬ 
age  repose  angle  was  53  deg.  A  third  series  of  tests  was  conducted 
using  pressure  fluctuations  to  simulate  the  turbulent  fluctuations  that 
occur  when  water  flows  over  the  riprap.  A  pressure  transducer  was 
installed  flush  with  the  sloping  side  on  which  the  riprap  was  placed. 
Measurements  of  pressure  were  taken  for  flow  conditions  close  to  the 
conditions  that  resulted  in  failure  of  the  riprap.  With  the  test  sec¬ 
tion  submerged  but  without  flow,  a  variable  speed  vibrator  was  attached 
to  the  flume  sidewall.  The  speed  of  the  vibrator  was  varied  until  the 
amplitude  of  the  measured  fluctuation  was  approximately  equal  to  the 
maximum  amplitude  measured  under  flowing  water  conditions.  This 
vibrator  speed  was  used  in  all  subsequent  angle  of  repose  tests.  The 
third  series  with  pressure  fluctuations  resulted  in  an  average  repose 
angle  of  53  deg.  The  vibrator  resulted  in  higher  frequency  fluctuations 
than  did  the  flowing  water  condition  but  the  amplitudes  were  similar. 

The  predictive  technique  of  Carter,  Carlson,  and  Lane  (1953)  was 
again  tested  against  the  experimental  data  using  the  measured  angle  of 
repose  of  53  deg.  Results  given  in  Figure  4.20  show  a  much  better  com¬ 
parison  between  predicted  and  observed  values  when  the  repose  angle  of 
53  deg  is  used  in  the  Carter,  Carlson,  and  Lane  equation,  Equation  2.16. 

Additional  tests  were  conducted  to  determine  why  the  measured 
repose  angle  was  significantly  higher  than  that  predicted  by  existing 
techniques  (Anderson,  PaintaL,  and  Davenport  1970).  These  tests  were 
conducted  to  determine  the  effects  of  revetment  height,  bank  smoothness, 
and  revetment  thickness.  Results  shown  in  Table  4.14  were  plotted  in 
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Figure  4.20.  Tractive  force  ratio  K  versus  0  (Data  from 

Table  4.13) 

Figure  4.21.  The  relative  height  of  the  revetment  is  defined  ns  the 
length  along  the  slope  L  divided  by  the  average  riprap  size  d^  . 
Revetment  height  Lg  was  included  to  determine  if  a  50-ft-high  channel 
bank  is  less  stable  or  has  a  different  angle  of  repose  than  a  10-ft-high 
channel  bank.  Also  shown  on  Figure  4.21a  is  the  repose  angle  for 
crushed  rock  from  Anderson,  Paintal,  and  Davenport  (1970).  These  results 
show  that  revetment  height  and  thickness  have  a  significant  effect  on 
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the  angle  of  repose.  Surface  smoothness  was  tested  by  comparing  the 
repose  angle  for  a  smooth  piece  of  marine  plywood  to  that  of  a  surface 
having  sand  glued  to  the  marine  plywood.  The  smooth  surface  yields  a 
slightly  higher  value  than  the  sand  surface.  The  difference  is  small 
and  surface  roughness  is  not  considered  to  be  a  big  factor  in  angle  of 
repose  for  the  two  surfaces  used  in  these  tests. 

The  California  Division  of  Highways  (1970)  uses  a  repose  angle  of 

70  deg  in  the  predictive  equations.  Blodgett  and  McConaughy  (1986) 

report  that  this  was  based  on  tests  in  which 

They  constructed  a  model  streambank  on  which  small  stones 
were  arranged  as  riprap,  and  underlying  stones  were  cemented 
in  a  plaster  of  paris  base.  The  side  slope  was  increased 
until  the  first  outer  stone  was  displaced.  It  was  determined 
that  65°  70°  was  the  maximum  angle  attained  before  a  stone 

fell  out. 

Miller  and  Byrne  (1965)  found  the  angle  of  repose  of  a  single  sand  grain 
on  a  fixed  rough  bed  to  be  as  high  as  70  deg  when  the  fixed  rough  bed 
particles  were  equal  in  size  to  the  single  sand  grain.  Both  the 
California  study  and  Miller  and  Byrne  show  that  surface  roughness 
becomes  important  when  the  underlying  material  size  becomes  large  rela¬ 
tive  to  the  size  of  the  riprap.  Hudson  (1958)  did  not  include  the  coef¬ 
ficient  of  friction  (angle  of  repose)  in  the  development  of  his  widely 
used  equations  for  the  design  of  quarrystone  cover  layers  subjected  to 
wave  attack.  He  cited  several  factors  that  presented  difficulty  in 
using  angle  of  repose.  Method  of  placement  was  one  of  the  factors  that 
caused  variation  in  the  repose  angle. 

The  following  results  summarize  angle  of  repose: 

1.  The  angle  of  repose  of  a  revetment  is  not  always  equal  to  the 
bulk  angle  of  repose  reported  in  the  literature. 


2.  The  angle  of  repose  of  a  revetment  is  affected  by  revetment 
height,  thickness,  method  of  placement,  and  possibly  other 
factors  that  were  not  investigated. 

Similar  to  Hudson  (1958),  this  investigation  will  omit  repose  angle  from 
the  analysis  of  side  slope  stability  and  incorporate  repose  angle 
effects  into  the  empirically  derived  coefficients. 

4.4.2  Velocity  Profiles  Over  Channel  Side  Slopes.  As  part  of  the 
CSU  Phase  IV  riprap  stability  tests,  velocities  were  measured  over  the 
IV: 2H  side  slope  in  a  straight  flume.  Results  from  tests  having  similar 
depth  were  averaged.  Velocities  were  made  dimensionless  by  dividing 
observed  point  velocity  by  the  average  velocity  of  a  single  vertical 
traverse  over  the  toe  of  the  slope,  and  depths  were  expressed  in  per¬ 
centage  of  the  total  depth.  Results  shown  in  Figure  4.22  indicate 
reduced  velocities  over  the  slope  and  that  the  influence  of  the  slope 
extends  out  from  the  toe  of  the  slope  approximately  0.5  times  the  depth 
of  flow.  The  measured  depth  at  the  toe  of  the  slope  was  generally 
95  percent  of  the  depth  in  the  horizontal  portion  of  the  channel.  This 
is  shown  in  Figure  4.22  where  the  cross  section  is  rounded  at  the  toe  of 
the  slope.  These  profiles  are  for  straight  channels  without  the  effects 
of  upstream  channel  curvature.  An  analysis  of  the  shear  distribution  of 
the  profiles  was  conducted  using  the  approach  given  in  Section  4.4.1. 

The  shear  stress  was  evaluated  relative  to  the  shear  stress  in  the  hori¬ 
zontal  portion  of  the  cross  section  at  X/D  =  -1.0  .  The  velocity  along 
the  channel  bottom  was  determined  at  a  distance  of  0.1D  above  the  bed, 


where  D  is  the  depth  at  X/D  =  -1.0  .  The  relative  shear  is 
determined  from 


T  at  g  =  -1.0  V  at  g  =  -1.0 


Results  from  the  three  profiles,  plotted  in  Figure  4.23,  show  that  the 
shear  stress  is  less  on  the  channel  side  slope  than  on  the  channel 
bottom. 

Another  series  of  velocity  measurements  was  conducted  in  a  curved 
channel  at  WES  to  determine  velocity  profiles  over  side  slopes  that  have 
strong  upstream  curvature  effects.  Profiles  were  measured  at  sta¬ 
tions  11.6,  16.6,  21.6,  65.0,  70.0,  and  75.0,  shown  in  Figure  4.24. 

These  stations  correspond  to  the  regions  of  maximum  velocity  over  the 
toe  of  the  slope.  Nondimensional  profiles  were  determined  (Figures  4.25 
and  4.26)  and  are  significantly  different  from  the  profiles  having  a 
straight  nnc-t-^eam  alignment.  These  curved  channel  profiles  show  a 
velocity  maximum  over  the  toe  of  the  slope,  with  the  maximum  located 
below  the  water  surface. 
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Figure  4.23.  Shear  stress  distribution,  1V:2H  side  slope,  straight 

channel 


Figure  4.24.  WES  curved  channel  model,  plan  view 

An  analysis  of  the  shear  in  the  straight  versus  curved  section  was 
conducted  using  the  same  analysis  used  in  Equation  4.24.  For  equal 
average  velocity  over  the  toe  of  the  slope,  the  maximum  stress  on  the 
curved  channel  side  slope  (located  at  X/D  =  0.5)  is  equal  to  approxi¬ 
mately  1.5  times  the  shear  stress  at  X/D  =0.5  in  the  straight  channel 
side  slope. 

4.4.3  Side  Slope  Stability  Tests.  Before  side  slope  stability 
tests  are  analyzed,  a  characteristic  velocity  and  depth  must  be 
selected.  This  velocity  and  depth  must  be  representative  of  the  condi¬ 
tions  on  the  side  slope  and  must  also  be  values  which  a  designer  has 
some  hope  of  determining  or  estimating.  Average  channel  velocity  is  the 
easiest  to  determine  but  not  very  representative  of  conditions  on  the 
side  slope.  Depth  and  average  velocity  over  the  toe  of  the  slope  will 
be  used  in  this  investigation  for  the  side  slope  stability  analysis. 
These  values  were  selected  based  on  the  two  requirements  stated  pre¬ 
viously  and  the  results  from  the  WES  curved  channel,  which  showed  that 


I 


the  maximum  velocity  in  the  cross  section  occurred  over  the  toe  of  the 
side  slope. 

Prior  to  the  CSU  Phase  IV  side  slope  stability  tests,  an  analysis 
was  conducted  of  the  Dorena  Dam  prototype  tests  reported  by  the  US  Army 
Engineer  District,  Portland  (1952).  These  tests  were  conducted  down¬ 
stream  of  Dorena  Dam  in  a  channel  having  a  grouted  riprap  bottom  and 
1V:2H  side  slopes  with  riprap  placed  to  a  thickness  of  ld^^  '  Results 
are  shown  in  Table  4.15  and  Figure  4.27.  The  curve  for  threshold  of 
incipient  failure,  thickness  of  ld^^  *  IV:  2H  side  slope,  and  straight 
channel  is 


based  on  the  Dorena  Dam  prototype  tests.  The  depth  and  velocity  over 
the  toe  of  the  slope  were  used  in  the  analysis.  The  velocities  were 
measured  for  several  tests,  and  these  were  used  to  estimate  the  velocity 
in  the  remaining  tests. 

In  the  CSU  Phase  IV  tests,  stability  was  determined  for  the 
following : 


d50 

Thickness 

in. 

dioo 

Table 

1.0 

1.33 

4.16 

1.0 

1.0 

4.17 

0.5 

1 .0 

cc 

Like  the  Dorena  Dam  tests,  these  tests  were  conducted  in  a  straight 
channel  without  upstream  curvature  effects. 

The  results  given  in  Tables  4.16-4.18  show  that  the  bottom  riprap 
fails  more  often  or  with  greater  severity  than  the  side  slope  riprap. 
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Figure  A. 27.  d^/D  versus  modified  Froude  number,  thickness  of 

1.0d.nn  ,  1V:2H  side  slope,  Dorena  Dam  prototype  tests  (Data 

from  Table  4.15) 

In  tests  29-35  (Table  4.16)  with  the  3-in.  d^^  riprap  placed  2  in. 

thick,  the  bed  was  stabilized  with  a  wire  screen  to  ensure  failure  on 
the  side  slopes.  Results  from  tests  with  a  thickness  of  ld^Q  are 
shown  in  Figure  4.28.  A  failure  point  from  test  21  is  located  to  the 
left  of  the  incipient  failure  line.  This  test  had  a  total  failed  area 
of  less  than  0.1  sq  ft.  Due  to  the  small  failure  area  and  the  position 
of  this  point  relative  to  several  stable  runs,  this  point  was  not  used 


in  the  determination  of  the  best-fit  line.  The  curve  for  threshold  of 
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Figure  4.28.  d^/D  versus  modified  Froude  number,  thickness 


of  ld,„„  ,  1V:2H  side  slope,  CSU  Phase  IV  (Data  from 
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Tables  4.17  and  4.18) 


incipient  failure,  a  thickness  of  ld10o  ’  1V:2H  side  slope’  and 


straight  channel  is 


So  =  4  /  yw  y/z  jl 
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based  on  the  CSU  Phase  IV  tests.  This  relation  is  in  close  agreement 


with  the  Dorena  Dam  prototype  test  results. 


Results  from  the  CSU  tests  with  a  thickness  of  1.33d1()0  are  shown 


in  Figure  4.29.  The  effects  of  thickness  for  side  slope  riprap  are 
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Figure  4.29.  d^/D  versus  modified  Froude  number,  a 


thickness  of  1.33d1Q0  ,  IV: 2H  side  slope,  CSU 


Phase  IV  (Data  from  Table  4.16) 


From  Figure  4.15  the  effect  of  thickness  for  bottom  riprap  is 


CU(N  =  1.33) 


0.225 


C15(N  =  1.00)  0.30 


=  0.75 


which  is  essentially  the  same  thickness  effect  as  the  side  slope  riprap. 


Mote  that  these  thickness  results  apply  only  to  gradations  having  a 


^85^15  rat^°  2-2.3. 


A  limited  series  of  stability  tests  was  conducted  in  the  WES  curved 


channel  facility  shown  in  Figure  4.24.  These  tests  were  limited  in  the 


sense  that  only  a  narrow  range  of  rock  size  could  be  failed  while  main¬ 


taining  a  high  enough  Reynolds  number  based  on  the  limitations  given  in 


section  3.5.  Results  are  shown  in  Table  4.19  and  Figure  4.30.  With  so 


few  data  points,  the  basic  relation  given  by  Equation  4.21  is  used  to 


define  the  slope  of  the  power  function.  The  relation  describing  the 
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Figure  4.30.  d^/D  versus  modified  Froude  number, 
thickness  of  »  1V:2H  side  slope,  WES  curved 

channel  (Data  from  Table  4.19) 

threshold  incipient  failure  for  a  thickness  of  ld100  »  IV: 2H  side  slope, 
and  curved  channel  is 


=  0.29 
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based  on  the  WES  curved  channel  model  using  depth  and  average  velocity 
over  the  toe  of  the  side  slope. 

This  model  derived  relationship  can  be  compared  to  the  prototype 
data  of  Blodgett  and  McConaughy  (1986).  These  data  were  taken  mainly 
from  curved  channels.  Since  side  slope  angle  is  generally  considered  to 
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be  significant,  only  sites  with  the  cotangent  of  the  side  slope  angle 
between  1.8  and  2.2  were  used  in  this  analysis.  To  best  estimate  the 
average  velocity  and  depth  over  the  toe  of  the  slope,  the  maximum  depth 
and  maximum  velocity  were  used  in  the  analysis.  In  cases  where  the 


maximum  velocity  was  not  measured,  the  relation 

V  =  1.53V  4.28 

max  avg 

was  derived  from  the  Blodgett  and  McConaughy  data  and  used  to  estimate 
maximum  velocity.  A  similar  approach  was  used  to  estimate  d^  for 
measurements  where  sufficient  data  were  not  given.  The  values  used  in 
the  analysis  are  shown  in  Table  4.20  and  Figure  4.31.  Only  sites  with 
particle  erosion  or  no  damage  were  considered  in  the  analysis.  Blanket 
thicknesses  are  not  given,  and  these  results  are  considered  applicable 
to  a  thickness  of  ld1Q0  since  most  prototype  sites  are  constructed  to 
this  thickness.  The  incipient  failure  curve  shown  in  Figure  4.31  is 
drawn  to  the  right  of  two  failure  points  (measurements  6  and  7) .  The 
site  of  these  measurements  was  a  channel  curved  only  18  deg.  The 
failure  for  measurement  6  was  on  the  upper  6  ft  of  the  channel  side 
slope,  which  is  unusual.  Blodgett  and  McConaughy  state  that  the 
velocity  for  measurement  6  "may  have  been  greater  than  estimated."  The 
data  for  measurement  7  show  a  relatively  low  velocity  but  an  extremely 
high  shear  stress.  The  high-water  profile  shows  some  unusual  conditions 
such  as  an  adverse  water-surface  slope  over  the  point  of  failure.  Con- 
slueLxug  these  problems  and  the  proximity  of  points  6  and  7  to  other 
stable  points,  these  failure  points  were  not  considered  in  the  analysis. 
The  resulting  threshold  of  incipient  failure  curve  for  IV: 2H  side  slope, 
a  thickness  of  ld^^  ,  and  curved  channels  is  the  same  as  that  proposed 
for  bottom  riprap  or 
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CHAPTER  5 
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VELOCITY.  ET/SEC 


a.  Velocity  versus  depth  versus 


for  unit  stone  weight  = 


165  pcf,  thickness  =  Id 


THICKNESS 


Thickness  adjustment  for 


d  /d  -  2. 0-2. 3 
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c.  Correction  for  unit 
stone  weight 


Figure  5.1.  Sizing  nomograph  for  riprap 


5.2  SIZING  NOMOGRAPH 


A  sizing  nomograph  of  Equation  5.1  is  shown  in  Figure  5.1.  In 


Figure  5.1a,  the  relationship  of  average  velocity  in  the  vertical, 


depth,  and  d  n  are  given  for  a  specific  weight  of  165  pcf  and  a 


blanket  thickness  of  ld^^  ’  ^&ure  5.1b,  the  adjustment  for 


thickness  is  given  for  gradations  having  dg^/d^  =  2.0-2. 3  ,  which  is 


similar  to  the  gradations  given  in  OCE  (1971).  In  Figure  5.1c,  the  cor 


rection  for  unit  weight  of  rock  is  given. 


5.3  DESIGN  APPLICATION 

This  design  procedure  is  based  on  the  premise  that  a  variety  of 
tools  are  available  for  estimating  the  average  velocity  in  the  vertical 
for  use  in  this  design  procedure.  The  average  velocity  in  the  vertical 
at  the  point  of  interest  is  used,  not  average  cross-section  values.  The 
available  tools  for  determining  velocities  include  the  following: 

1.  Numerical  Models:  one-dimensional  water-surface  profile 
programs  and  multidimensional  models 

2.  Physical  models 

3.  Prototype  measurements 

4.  Analytical  techniques  such  as  the  California  Division  of 
Highways  (1970)  equation 

V  =  4/3 (Average  Channel  Velocity)  5.2 

which  should  be  limited  to  prismatic  channels,  and  the  analysis 
of  the  data  by  Blodgett  and  McConaugby  (1986)  which  gave  the 
relation 

V  =  1*^3 (Average  Channel  Velocity)  (4.28  bis) 

which  would  be  applicable  to  natural  channels. 

Analytical  techniques  that  need  to  be  developed  include 

V  =  /(average  channel  velocitv,  bend  radius/ 
max 

water-surface  width,  channel  shape  (natural 
or  prismatic),  side  slope  angle,  aspect  ratio, 
bend  angle,  different  bed  and  bank  roughness) 

5.4  EXAMPLE  DESIGN 

Determine  size  of  side  slope  riprap  for  the  design  problem  at 
pinole  Creek  given  in  Blodgett  and  McConaughy  0986)  having  the  fol¬ 


lowing  conditions: 


Average  channel  velocity  =  7.7  ft/sec 
Average  depth  =  4.8  ft 
Maximum  depth  =  7.7  ft 
Yg  =  175  pcf 

Curved  channel  (radius/width  =  2.5) 

Water-surface  width  =  61  ft 

Thickness  =  ld^gg 

Cotangent  of  side  slope  angle  =  2 

As  in  most  riprap  design  problems,  only  the  average  channel  velocity  is 
known.  The  maximum  average  velocity  in  the  vertical  over  the  toe  of  the 
outer  bank  can  be  estimated  by 

V  =  1.53V  =  11.78  ft/sec  (4.28  bis) 

max  avg 

This  velocity  and  the  maximum  depth  are  used  in  equation  5.1.  The  size 
required  for  stability  is  d^g  =  0.64  ft  .  At  Pinole  Creek  prototype, 
riprap  having  a  d^g  of  0.45  ft  and  a  unit  stone  weight  of  178  pcf 
failed  under  the  given  hydraulic  conditions. 

Using  OCE  (1971)  gradations  given  in  Table  4.21,  a  blanket  thick¬ 
ness  of  18  in.  provides  a  d^g  (minimum)  of  0.73  ft  for  a  unit  stone 
weight  of  175  pcf.  A  blanket  thickness  of  15  in.  cannot  be  used  because 
the  d^g  (minimum)  of  0.61  ft  is  less  than  0.64  ft. 

For  comparison,  OCE  (1970  and  1971)  riprap  sizing  guidance  using  a 
constant  Shields  coefficient  and  the  logarithmic  velocity  relations 
results  in  a  d^g  of  1.17  ft.  The  24-in.  blanket  thickness  given  in 
Table  4.21  for  a  unit  stone  weight  of  175  pcf  provides  a  d^g  (minimum) 


5.5  SUMMARY  OF  LIMITATIONS 


This  design  procedure  is  limited  to  the  following  conditions: 

1.  Straight  and  curved  open  channels  that  are  not  immediately 
downstream  of  a  structure  that  creates  a  hydraulic  jump. 


2.  Channel  bottoms  and  channel 
IV: 2H. 

3.  Slopes  less  than  2  percent, 

4.  Froude  number  less  than  1.2. 

5.  Ratio  of  flow  depth  to  d^^ 

6.  For  thickness  equal  ld^^  * 
4.6.  For  thickness  greater 

7.  Angular  rock. 


side  slopes  less  than  or  equal  to 

no  overtopping  embankment  flows. 

riprap  size  from  4  to  50. 
dgs /d i 5  ,  less  than  or  equal  to 
than  *  d85/,d15  ^rom  2-°~2.3. 


CHAPTER  VI 


CONCLUSIONS  AND  RECOMMENDATIONS 


This  investigation  has  shown  that  a  constant  Shields  coefficient 
and  the  logarithmic  velocity  laws  are  not  applicable  to  such  high  rela¬ 
tive  roughness  problems  as  riprap  design. 

The  critical  velocity  relation  developed  in  this  study  for  the 
threshold  of  incipient  failure  of  riprap  is 


^30  _  /  Yw  \1/2 
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This  relation  was  developed  for  straight  channel  bottoms,  straight  chan¬ 
nel  side  slopes  of  1V:2H,  and  curved  channel  side  slopes  of  1V:2H. 
Average  velocity  in  the  vertical  at  the  point  of  interest  is  used,  not 
average  cress-sectional  values.  A  relation  of  this  form  was  first  pro¬ 
posed  by  Neill  (1967)  and  Bogardi  (1968). 

This  critical  velocity  relation  was  compared  to  the  two  most  common 

riprap  sizing  methods:  (1)  critical  shear  stress  using  a  constant 

2 

Shields  coefficient  and  (2)  Isbash  type  relations  (d,_Q  =  C^V  ).  This 
critical  velocity  relation  has  an  exponent  that  falls  halfway  between 
these  two  methods. 

Riprap  gradation  uniformity  was  shown  to  affect  riprap  stability  if 
d.(j  ’  s  used  in  the  analysis.  Use  of  particle  size  d^g  in  the 
stability  relations  eliminates  the  effects  of  gradation  uniformity  for 
rip  at-  thickness  of  ld^  g  and  is  used  as  the  characteristic  size  in 
t  L  :  -  i  nves  t  i  g.i  t  i  on  . 
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Riprap  thickness  was  shown  to  have  a  significant  affect  on  riprap 

stability  for  riprap  gradation  having  doc/d,r  of  2. 1-2.3. 

od  io 

Within  the  range  tested,  riprap  shape  did  not  have  a  significant 
effect  on  riprap  stability.  Gradations  having  18  percent  elongated 
particles  (2,/ b  >  3)  exhibited  the  same  stability  as  gradations  not 
having  elongated  particles. 

Existing  side  slope  relations  used  in  the  critical  shear  stress 
equation  overestimate  the  decrease  in  stability  that  occurs  when  a 
particle  is  placed  on  a  sloping  bank.  This  was  demonstrated  in  two 
ways : 

1.  Comparison  with  the  Hudson  (1958)  wave  equation  that  showed  the 
effects  of  side  slope  angle  are  more  significant  in  channel 
flow  than  in  wave  attack. 

2.  Stability  tests  on  sloping  sides  conducted  in  this  investiga¬ 
tion.  The  existing  side  slope  stability  relations  matched  the 
observed  data  when  a  repose  angle  of  53  deg  was  used  in  the 
analysis  instead  of  the  commonly  used  40  deg.  This  led  to  a 
series  of  repose  angle  tests  which  suggested  that  repose  angle 
varies  with  revetment  height,  riprap  thickness,  surface  tex¬ 
ture,  and  placement  method.  Revetment  height  is  important 
because  the  higher  the  bank,  the  greater  the  amount  of  material 
being  supported  by  the  rock  at  the  toe  of  the  slope.  Because 
of  these  difficulties,  repose  angle  was  not  used  in  the  criti¬ 
cal  velocity  relation  and  was  included  in  the  empirical  coef¬ 
ficients  just  as  Hudson  did  in  his  wave  equation. 

Comparison  of  velocity  profiles  over  channel  side  slopes  in 
straight  and  curved  reaches  shows  that  for  the  same  average  velocity 


over  the  toe  of  the  side  slope,  the  velocity  and  shear  sticoo  on  the 
side  slope  are  significantly  higher  on  the  outer  bank  of  the  curved 
channel.  Side  slope  stability  tests  in  straight  channels  cannot  be  used 
in  channel  bends.  The  representative  velocity  used  herein  for  side 
slope  riprap  stability  is  the  average  velocity  in  the  vertical  over  the 
toe  of  the  side  slope.  Maximum  velocities  in  the  curved  channel  tests 
occurred  over  the  toe  of  the  side  slope. 

Using  depth  and  velocity  over  the  toe  of  the  side  slope,  the 
threshold  of  incipient  failure  of  IV: 2H  side  slope  riprap  in  straight 
channels  is  described  by 
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based  on  model  and  prototype  data  for  thickness  =  ldJOg. 

For  IV: 2H  side  slope  riprap  in  curved  channels 
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based  on  model  and  prototype  data  for  thickness  of  ld100  •  This  rela¬ 
tion  was  also  found  applicable  to  bottom  riprap  in  straight  channels. 

Since  these  relations  define  the  threshold  of  incipient  failure, 
safety  factors  must  be  determined  before  they  can  be  used.  A  common 
problem  in  the  design  of  riprap  is  the  addition  of  safety  factors  at 
each  step  in  the  design  procedure.  A  safety  factor  of  1.2  times  the  d^Q 
riprap  size  given  by  the  threshold  of  incipient  failure  curve  is  used  in 
the  sizing  nomograph  (Figure  5.1)  developed  in  this  investigation.  The 
designer  can  easily  use  other  safety  factors  and  apply  them  to  the 
incipient  failure  relations. 


During  this  study,  the  following  areas  were  identified  as  needing 


additional  study  relative  to  open  channel  riprap  design: 

1.  Effects  of  blanket  thickness  for  gradations  other  than  those 
studied  in  this  investigation. 

2.  Effects  of  riprap  shape  outside  the  range  covered  in  this 
study,  including  the  effects  of  surface  texture  such  as 
stability  of  cobble  particles. 

3.  Side  slope  stability  tests  of  1V:1.5H  and  1V:3H. 

4.  Determining  repose  angle  of  riprap  revetment  so  that  it  can  be 
included  in  the  design  procedure. 

5.  Effect  of  revetment  side  slope  height  on  stability.  Side  slope 
riprap  in  shallow  channels  may  be  much  more  stable  than  in  deep 
channels  due  to  the  amount  of  material  being  supported  by  the 
toe  of  the  slope. 

6.  Analytical  methods  for  determining  velocity  in  straight  and 
curved  channels  for  use  in  riprap  sizing. 

7.  Using  the  experience  of  others  involved  in  riprap  design  to 
better  define  appropriate  safety  factors. 
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0.025 

0.029 

3/8-774 

0.00600 

3.03 

0.377 

0.019 

0.025 

0.029 

3/8-774 

0.00400 

2.88 

0.487 

0.019 

0.025 

0.029 

3/8 -7<4 

0.00500 

3.00 

0.461 

0.019 

0.025 

0.029 

3/8-774 

( Cont inued) 


Velocity  based  on  discharge/area, 

S  =  stable;  F  =  failed;  ?  =  unknown. 

^  Stable  but  tested  for  a  short  duration  compared  to  the  other 
Width/depth  <  5. 


tests 


4WJ. 


.  '  «  M  -  *  * 


Table  4.1  (Concluded') 


Energy  Average  Average 
Slope  Velocity  Depth 

c  .  I  r  .  r- .  /  ^  j"  . 


ft/ft  ft/sec 


d30/ft  d50/ft  d90/ft  Size,  in. 


Stable  or 
Failed  or 

Unknown*5 


Table  4.2 


CSU  Phase  I, 

Test  Results  for  dg^/d^  =  2.8,  Thickness  =  2d,.Q  =  Id^Q, 
d50  =  1.87  in.,  Yg  =  170  pcf,  Angular  Particles,  Does  Not 


Meet  Corps  Shape  Guidance 


Orifice 

Average 

Average 

Run 

Discharge 

FI  "roe 

Velocity 

Depth 

Stable  (S)  or 

No. 

cf  s 

Slope 

fps 

ft 

Failed  (F) 

1 

25 

0.00852 

3.62 

0.825 

S 

2 

0.01378 

4.34 

0.687 

S 

3 

0.01667 

4.62 

0.653 

F 

4 

0.01973 

5.11 

0.629 

S 

5 

0.00761 

4.78 

1.218 

S 

6 

0.01089 

5.24 

1.110 

S 

7 

0.01451 

5.92 

1.032 

F 

8 

0.01766 

6.01 

1.051 

F 

9 

0.01089 

5.35 

1.151 

F 

10 

75 

0.00537 

5.39 

1.8683 

S 

11 

0.00769 

5.38 

1 .71 6a 

S 

12 

0.01025 

5.92 

1.509 

F 

13 

0.00894 

5.64 

1.6043 

F 

14 

0.00769 

5.39 

1.7003 

S 

15 

100 

0.00420 

5.39 

2.325a 

S 

16 

0.00601 

5.32 

2. Ill3 

S 

17 

0.00801 

F 

18 

0.00699 

5.86 

2.012° 

F 

19 

0.00601 

5.68 

2.0623 

S 

a 


Width/depth  <  5  . 
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de.  =  0.5  in.,  v  =  166  pcf 
50  s 


Nominal 
Discharge 
cf  s 


Flume 

Slope 


Average 

Velocity 


Average 

Depth 

ft 


Stable  (S)  or 
Failed  (F) 


a 


Width/depth  <  5 


Table  4.4 


CSU  Phase  III, 

Test  Results  for  dg^/d^  =  2.1,  Thickness  =  2d^  =  1.4d^QQ, 
d^Q  =  1  in.,  ys  =  167  pcf.  Shape  Characteristics 


Meeting  Corps  Guidance 


Nominal 

Average 

Average 

Run 

Discharge 

Flume 

Velocity 

Depth 

Stable  (S)  o 

No. 

cf  s 

Slope 

fps 

ft 

Failed  (F) 

33 

25 

0.00998 

4.40 

0.684 

S 

34 

25 

0.01088 

4.51 

0.703 

S 

35 

25 

0.01186 

4.72 

0.672 

S 

36 

25 

0.01337 

4.95 

0.618 

F 

37 

25 

0.01204 

4.77 

0.651 

F 

32 

50 

0.00558 

4.90 

F 

41 

50 

0.00475 

4.71 

1.353 

F 

38 

75 

0.00402 

5.02 

1.832a 

F 

39 

75 

0.00377 

1 . 842a 

F 

40 

75 

0.00345 

4.84 

1.91 8a 

S 

42 

100 

0.00314 

4.97 

2 . 3  7 1 a 

S 

43 

100 

0.00403 

4.90 

2.4153 

S 

44 

100 

0.00436 

5.20 

2.210a 

F 

45 

0.00354 

5.09 

2. 332a 

F 

3  Width/depth  <  5  . 


Table  4.5 
CSU  Phase  III, 

Test  Results  for  =  2.3,  Thickness  =  2d^  =  1.4djQQ, 

d^Q  =  2  in.,  yg  =  165  pcf.  Shape  Characteristics 


Meeting  Corps  Guidance 


Table  4.6 


Width/depth  <  5 


Table  4.7 


CSU  Phase  III, 


Test  Results  for  dQC/d.c  =  2.3,  Thickness  =  3dcn  =  2.1d,__, 

oj  I  j  !>U  1UU 


d^  =  2  in.,  =  165  pcf.  Shape  Characteristics 


8 


Width/depth  <  5 


Meeting  Corps  Guidance 


Nominal 

Average 

Averagi 

Run 

Discharge 

Flume 

Velocity 

Depth 

No. 

cfs 

Slope 

fps 

ft 

79 

25 

0.01180 

4.42 

0.710 

80 

25 

0.01870 

5.17 

0.607 

81 

50 

0.01205 

5.90 

1.068 

82 

50 

0.01544 

6.47 

0.966 

83 

50 

0.01724 

6.76 

0.928 

84 

50 

0.01879 

6.61 

0.970 

85 

75 

0.00898 

6.19 

1.519 

86 

75 

0.01095 

6.58 

1.414 

87 

75 

0.01206 

6.63 

1.423 

88 

75 

0.01359 

6.88 

1.372 

89 

75 

0.01565 

6.84 

1.399 

90 

100 

0.00866 

6.97 

1 . 808a 

91 

100 

0.00938 

6.96 

1 . 796a 

92 

100 

0.01084 

7.39 

1 . 7 1 1 a 

93 

100 

0.01189 

7.44 

1 . 698a 

94 

100 

0.01300 

8.02 

1.572 

Stable  (S)  or 
Failed  (F) 


O  ^  CO 


Table  4.8 


CSU  Phase  I, 

Test  Results  for  dg^/d^  =  3.9,  Thickness  =  2d^  =  ld^^. 


d_,.  =  3  in.,  y  =  170  pcf.  Angular  Particles,  Does  Not 
jU  s 


Meet  Corps  Shape  Guidance 


Orifice 

Average 

Average 

Run 

Discharge 

Flume 

Velocity 

Depth 

Stable  (S)  or 

No. 

cf  s 

Slope 

fps 

ft 

Failed  (F) 

1 

25 

0.02000 

4.51 

0.655 

S 

2 

50 

0.01544 

5.55 

1.064 

S 

3 

0.02000 

1.026 

S 

4 

75 

0.01500 

1.348 

F 

5 

0.01719 

6.55 

1.363 

F 

6 

0.01500 

6.72 

1.387 

F 

7 

0.01291 

6.41 

1.401 

S 

8 

100 

0.01009 

6.14 

1.831a 

s 

9 

0.01343 

6.54 

1 . 703a 

F 

10 

0.01172 

6.37 

1 . 825a 

F 

i 

i 


a 


Width/depth  <  5  . 


d50  =  1  ln*’ 


Yg  =  166  pcf 


Width/depth  <  5 


sws 


was 


v>- 

V,  ^./N  " 


Not  Meeting  Corps  Guidance 


Nominal 
Discharge 
cf  s 


Flume 

Slope 


Average 

Velocity 

fps 


Average 

Depth 

ft 


Stable  (S)  or 
Failed  (F) 


Width/depth  <  5 


Sws? 


vr  * 


\  ^  “V«L^ 


Table  4.11 

WES  Tilting  Flume,  Phase  II 


Bottom 

Slope 

ft/ft 


Average 

Velocity 

ft/sec 


Average 
Depth,  ft 


d50/ft 


d90/£t 


NOTE:  These  resistance  tests  were  conducted  using  riprap  with  one-third 

(by  weight)  3/4-1/2  in.,  one-third  1/2-3/8  in.,  and  one-third  3/8  in.-//4. 
The  stone  had  thickness  of  ld.^  .  Velocity  was  based  on  the  average 
of  two  vertical  profiles  taken  at  the  flume  center  line.  Stability  was 
not  studied  in  these  tests. 

3  Width/depth  <  5  . 
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Table  4.12 

WES  Tilting  Flume,  Phase  III 


Bottom 

Slope 

ft/ft 

Average 
Velocity , 
ft/ sec 

Average 
Depth,  ft 

d50/fta 

d90^ft 

0.00300 

1.94 

0.346 

0.036 

0.055 

0.00400 

2.09 

0.320 

0.036 

0.055 

0.00500 

2.26 

0.296 

0.036 

0.055 

0.00600 

2.38 

0.279 

0.036 

0.055 

0.00700 

2.99 

0.407 

0.036 

0.055 

0.01600 

2.72 

0.207 

0.036 

0.055 

0.00200 

2.10 

0.582 

0.036 

0.055 

0.00250 

2.37 

0.538 

0.036 

0.055 

0.00300 

2.48 

0.511 

0.036 

0.055 

0.00350 

2.60 

0.490 

0.036 

0.055 

0.00400 

2.69 

0.466 

0.036 

0.055 

0.00200 

1.73 

0.423 

0.036 

0.055 

0.00300 

1.96 

0.369 

0.036 

0.055 

0.00400 

2.16 

0.338 

0.036 

0.055 

0.00500 

2.28 

0.314 

0.036 

0.055 

0.00600 

2.40 

0.298 

0.036 

0.055 

0.00200 

2.06 

0.624 

0.036 

0.055 

0.00250 

2.16 

0.585 

0.036 

0.055 

0.00300 

2.33 

0.548 

0.036 

0.055 

0.00350 

2.54 

0.513 

0.036 

0.055 

0.00300 

1.81 

0.353 

0.051 

0.076 

0.00400 

2.07 

0.323 

0.051 

0.076 

0.00500 

2.26 

0.298 

0.051 

0.076 

0.00600 

2.42 

0.287 

0.051 

0.076 

0.00200 

2.12 

0.580 

0.051 

0.076 

0.00300 

2.30 

0.531 

0.051 

0.076 

0.0040C 

2.55 

0.488 

0.051 

0.076 

0.00500 

2.76 

0.446 

0.051 

0.076 

0.01300 

3.50 

0.379 

0.051 

0.076 

0.01200 

3.70 

0.425 

0.051 

0.076 

NOTE:  These  resistance  tests 

were  conducted  with  stone  having 

a  thick- 

ness  equal  to  Id 


Velocity  was  based  on  the  average  of  four  verti- 


W-WJVJtiJ.  WW  1  u  1  An  *  *  tJ.WWJ.UJ  «UvJ  wa  J  WVJ  Wll  tuv  W  V  V.4.  VI.  A.  VU1- 

cal  velocity  profiles.  Stability  was  not  studied  in  these  tests. 


Gradation  used  for 
Gradation  used  for 


0.036  ft  was  same  as  WES  Phase  II  tests. 
d,.Q~  =  0.051  ft  was  one-third  3/4-1  in., 


50 


one-third  1/2-3/4  in.,  and  one-third  3/8-1/2  in. 
Width/depth  <  5  . 


'  ^  ..A  *  At  K.  -  A  ft&.JL  M 


-% 


a 


Table  4.13 

Side  Slope  Versus  Critical  Bottom  Velocity 


Critical  Bottom 

Slope 

Velocity,  ft/sec 

Horizontal 

2.64 

2.58 

2.53 

IV :  4H 

2.63 

2.58 

2.63 

IV: 2 . 75H 

2.51 

2.58 

2.61 

2.61 

IV:  2H 

2.41 

2.44 

2.41 

2.46 

IV: 1.5H 

2.15 

2.20 

IV: 1 .25H 

2.06 

2.06 

1.87 

1.94 

1.91 

2.00 

*  *  *  V 


iXaVaa  al  A  hx  m  aA  A  aVAAA 


>V  .\  ‘»  A  V  V  .’ 

^  ^  -**  V  a^v\a.'-1.'.a\a 


Pressure 

Fluctuations  Submerged 


Revetment 

Thickness 

rtA 


Revetment 

Height 

L„ 


Average 
Angle  of 
No.  of  Repose 
Surface  Tests  de 


l!* 


Ht'll, 


Lt.J 


rJT’>'w>-T. 


wuvw. 


J"  ~s~  -r- 
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Table  4.15 


Analysis  of  Dorena  Dam  Prototype  Data 


IV; 2H  Side  Slope,  Straight  Channel 


I 


Sta- 

Right 

or 

Left 

Grada- 

Esti¬ 

mated 

V  Over 

Depth 

Over 

(y  Y/2 

W  1  JL  d30 

Failed 

or 

1  V  -Y  1  Vpfi  JU 

Test  tion 

Bank 

tion 

Toe 

Toe 

\Ys  Yw/  V8°  D 

Stable 

3 

2+80 

L 

A 

7.1 

7.3 

0.34 

0.067 

S 

3 

2+80 

R 

B 

7.1 

7.3 

0.34 

0.084 

S 

3 

3+46 

L 

A 

11.5 

6.6 

0.58 

0.074 

S 

3 

3+46 

R 

B 

10.5 

6.6 

0.53 

0.092 

S 

3 

3+90 

L 

C 

14.5 

6.0 

0.77 

0.095 

F 

3 

3+90 

R 

D 

14.5 

6.0 

0.77 

0.120 

S 

3 

4+12 

L 

C 

13.0 

6.7 

0.65 

0.085 

S 

3 

4+12 

R 

D 

13.0 

6.7 

0.65 

0.107 

S 

3 

4+75 

L&R 

D 

15.0 

6.3 

0.77 

0.114 

F 

4 

3+46 

L 

A 

10.9 

6.6 

0.55 

0.074 

S 

4 

3+46 

R 

B 

10.9 

6.6 

0.55 

0.092 

S 

4 

4+12 

L 

C 

13.0 

6.7 

0.65 

0.085 

S 

4 

4412 

R 

D 

13.0 

6.7 

0.65 

0.107 

S 

Rock  Characteristics: 


Gradation 

d30 

ft 

d100 

ft 

Thickness 

dioo 

A 

0.49 

1.08 

1.08 

B 

0.61 

1.28 

1.02 

C 

0.57 

1 

? 

D 

0.72 

1.79 

1.12 

W  C/3  fa 


Test  Results  for  doc/d.,  =  2.3,  Thickness  =  1.3d.,.,,,  dc„ 
oo  10  100  50 


Y  =  167  pcf,  Shape  Characteristics  Not  Meeting 


Corps  Criteria 


Run 
No. 3 

Discharge , 
cf  s 

Water- 
Surface 
Slope , 
ft/ft 

Average 
Velocity 
Over  Toe  of 
Slope 
ft/sec 

Average 

Depth, 

Over 

Toe , 
ft 

Stable  (S)  or 
Failed  (F) 

Bottom  Side  Slope 

1 

15 

0.00768 

3.44 

0.81 

S  S 

3 

15 

0.00929 

3.59 

0.76 

S  S 

4 

15 

0.01127 

3.73 

0.70 

F  S 

5 

15 

0.01077 

4.15 

0.64 

F  S 

6 

15 

0.00907 

4.18 

0.64 

F  S 

7 

15 

0.00957 

4.40 

0.62 

S  S 

12 

20 

0.00491 

3.91 

0.81 

S  S 

13 

20 

0.00804 

4.19 

0.78 

S  S 

14 

20 

0.00945 

4.41 

0.77 

S  S 

15 

20 

0.01074 

4.61 

0.77 

F  S 

16 

30 

0.00685 

4.86 

1.11 

S  S 

17 

30 

0.00723 

4.97 

1.09 

S  S 

18 

30 

0.00677 

5.13 

1.07 

S  s 

19 

30 

0.00796 

5.55 

1.00 

F  S 

20 

40 

0.00515 

4.87 

1.51 

S  S 

21 

40 

0.00595 

4.98 

1.42 

S  S 

22 

40 

0.00536 

5.25 

1.36 

S  S 

23 

40 

0.00669 

5.61 

1.28 

S  F 

25 

40 

0.00729 

5.76 

1.26 

F  F 

8 

50 

0.00547 

5.47 

1.57 

S  F 

9 

50 

0.00498 

5.46 

1.58 

S  F 

10 

50 

0.00526 

5.45 

1.60 

F  S 

11 

50 

0.00292 

5.56 

1.64 

S  S 

(Continued) 


a 


Test  numbers  omitted  did  not  have  velocities  measured  over  toe 


Hjhrj  ^  (/)(/)  (/)  MMW  W  W  W  W  W 


Table  4.16  (Concluded) 


Si 


scharge , 
cf  s 

Water- 
Surface 
Slope , 
ft/ft 

Average 
Velocity 
Over  Toe  of 
Slope 
ft/sec 

Average 

Depth, 

Over 

Toe , 
ft 

Stable 

Failed 

Bottom 

(S)  or 
(F) 

Side  Slope 

50 

0.00451 

5.77 

1.61 

— b 

S 

50 

0.00449 

5.89 

1.50 

— b 

F 

40 

0.00560 

5.71 

1.28 

— b 

F 

20 

0.00921 

4.15 

0.81 

— b 

S 

20 

0.01112 

4.64 

0.76 

— b 

S 

20 

0.01310 

4.77 

0.74 

— b 

F 

Test  numbers  omitted  did  not  have  velocities  measured  over  toe. 
Bottom  fixed  with  wire  mesh  to  ensure  side  slope  failure. 


V.v.v.v. 


>  V '."  WW  V  V.W  '/.V//  V 


I 


I 
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Table  4.17 
CSU  Phase  IV 

Test  Results  for  dg^/d^  =  2.3,  Thickness  =  Id^^,  d^^  =  1  in., 
Yg  =  167  pcf.  Shape  Characteristics  Not  Meeting 
Corps  Criteria 


Run 

No. 

Discharge , 
cf  s 

Water- 
Surface 
Slope , 
ft/ft 

Average 
Velocity 
Over  Toe  of 
Slope 
ft/sec 

Average 

Depth 

Over 

Toe , 
ft 

Stable 

Failed 

Bottom 

(S)  or 
(F) 

Side  Slope 

36 

15 

— b 

3.54 

0.76 

S 

S 

37 

15 

0.01002 

3.73 

0.68 

S 

S 

38 

15 

0.01165 

3.91 

0.64 

F 

S 

39 

15 

0.01090 

4.03 

0.69 

S 

S 

40 

20 

0.00900 

4.42 

0.80 

F 

S 

41 

20 

0.00832 

4.26 

0.84 

S 

S 

42 

30 

0.00711 

4.97 

1.09 

F 

F 

43 

30 

0.00482 

4.56 

1.20 

S 

S 

44 

30 

0.00649 

4.68 

1.16 

F 

FC 

45 

40 

0.00464 

4.93 

1.45 

F 

S 

46 

40 

0.00408 

4.62 

1.54 

S 

s 

47 

50 

0.00287 

4.93 

1.77 

s 

S 

48 

50 

0.00434 

5.27 

1.63 

F 

FC 

49 

50 

— b 

5.36 

1.63 

S 

F 

50 

50 

— b 

5 . 46 

1.61 

F 

F 

a 

b 

c 


Estimated  from  results  given  in  Table  4.16. 
Not  determined.  „ 


Failed  area  less  than  0.1  ft 
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Table  4.18 
CSU  Phase  IV 

Test  Results  for  dg^/d^  =  2.0,  Thickness  =  ld^^,  d.^  =  0.5  in 


167  pcf.  Shape  Characteristics  Not  Meeting 


Average 

Average 

Water- 

Velocity 

Depth 

Surface 

Over  Toe  of 

Over 

Stable  (S)  or 

Discharge , 

Slope , 

Slope 

Toe , 

Failed  (F) 

cfs 

ft/ft 

ft/sec 

ft 

Bottom  Side  Slope 

0,00203 

0.00269 

0.00207 

0.00197 

0.00375 

0.00295 

0.00400 

0.00347 

0.00242 

0.00234 

0.00221 

0.00206 

0.00322 

0.00270 

0.00240 

0.00250 

0.00241 

0.00170 

0.00268 

0.00159 

0.00158 


02 

S 

S 

93 

s 

S 

91 

F 

S 

Failed  area  less  than  0.10  ft 


-.’  V.-  ■*..  s_-  v-  O  O  ..- 


..  *  K  "  » 


Table  4.19 

Test  Results  From  WES  Curved  Channel  Model 


1 

Average  Velocity 

Average  Depth 

Number 

H 

Over  Toe  of 

Over 

of  6- 

H  Discharge 

Slope 

Toe , 

Stable  (S)  or 

Hour 

K.  cfs 

f t/seca 

ft 

Failed  (F) 

Runs 

For  stable  runs  this  was  the  maximum  average  velocity  in  the  vertical 
over  the  toe.  For  failure,  run  velocity  was  measured  at  the  location 
of  the  failure.  Failure  points  and  maximum  velocities  were  always 
between  stations  70  and  75. 


K*»»5 


v.v 

>,'•* 

v.v 

v.-y-y 


(From  Blodgett  and  McConaughy  (1986)) 


Measurement 

Number 


D 

max 


K  y/2 

v„«,  i— 

ft/sec  \  S  W  / 


V 


Riprap 

Thickness,  in 


_ Limits  of  Stone  Weight,  lb _ 

Percent  Lighter  by  Weight 
100  50  15 

Max  Min  Max  Min  Max  Min 


Stone  weight  limit  data  from  OCE  1971. 


Riprap 

Thickness,  in 


APPENDIX  A 


b 

bis 

B 


C1’C2,C3 
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d50,etc. 


d85/dI5 
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NOTATION 

effective  area  of  particle 
stone  breadth  or  thickness 
equation  is  repeated 
flume  bottom  width 
generic  coefficients 
Shields  coefficient 
drag  coefficient 
lift  coefficient 
flow  depth 

particle  size  of  which  a  certain  percent  is  finer  by 
weight 

gradation  uniformity 
lift  force 

shape  factor  and  surface  texture 
side  slope  factor 
universal  gravitational  constant 
wave  height 
tractive  force  ratio 
generic  coefficients 
stability  coefficient 
equivalent  sand  grain  roughness 
stone  length 
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Lg  =  length  along  channel  side  slope 

M,L,T  =  fundamental  dimensions  of  mass,  length,  and  time, 
respectively 

n  =  Manning's  roughness  coefficient 
N  =  blanket  thickness/d^Q 
R  =  hydraulic  radius 
S  =  energy  slope;  channel  slope 
=  shear  velocity  =  ^gDS  T/p 
V  =  average  flow  velocity 

=  cross-section  average  channel  velocity 
=  bottom  velocity 

=  critical  velocity  for  particle  on  horizontal  bed 
=  critical  velocity  for  particle  on  side  slopes 
=  local  velocity  at  distance  y 
W  =  unsubmerged  stone  weight 
W  =  submerged  stone  weight 

y  =  distance  above  origin  of  logarithmic  velocity  profile 

y  =  distance  below  top  of  roughness  element  to  origin  of 
°  profile 

a  =  bottom  angle  with  horizontal  in  flow  direction 

B  =  angle  of  inclination  of  drag  force  as  a  result  of 
secondary  motion 

=  specific  weight  of  stone 

Yw  =  specific  weight  of  water 

0  =  angle  of  side  slope  with  horizontal 

x  =  von  Karman  coefficient 


p  =  absolute  viscosity 
v  =  kinematic  viscosity 


fluid  density 
stone  density 

tractive  force  imposed  by  flowing  water 

critical  tractive  force  for  given  particle  size  on 
bottom 

critical  tractive  force  for  particle  on  side  slope 
angle  of  repose 
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